SU(2) gauge theory of the Hubbard model and application to the honeycomb lattice 
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Motivated by recent experiments on the triangular lattice Mott-Hubbard system k-(BEDT- 
TTF)2Cu2(CN)3, we develop a general formalism to investigate quantum spin liquid insulators 
adjacent to the Mott transition in Hubbard models. This formalism, dubbed the SU(2) slave-rotor 
formulation, is an extension of the SU(2) gauge theory of the Heisenberg model to the case of the 
Hubbard model. Furthermore, we propose the honeycomb lattice Hubbard model (at half-filling) as 
a candidate for a spin liquid ground state near the Mott transition; this is an appealing possibility, 
as this model can be studied via quantum Monte Carlo simulation without a sign problem. The 
pseudospin symmetry of Hubbard models on bipartite lattices turns out to play a crucial role in our 
analysis, and we develop our formalism primarily for the case of a bipartite lattice. We also sketch 
its development for a general Hubbard model. We develop a mean-field theory to describe spin 
liquids and some competing states, and apply it to the honeycomb lattice. On the insulating side 
of the Mott transition, we find an SU(2) algebraic spin liquid (ASL), described by gapless S = 1/2 
Dirac fermions (spinons) coupled to a fluctuating SU(2) gauge field. This result contrasts with that 
obtained via a U(l) slave-rotor approach, which instead found a U(l) ASL. That formulation does 
not respect the pseudospin symmetry, and is therefore not correct on the honeycomb lattice. We 
construct a low-energy effective theory describing the ASL phase, the conducting semimetal phase 
and the Mott transition between them. This physics can be detected in numerical simulations via 
the simultaneous presence of substantial antiferromagnetic and valence-bond solid correlations. In 
the SU(2) ASL, these observables have slowly-decaying fluctuations in space and time, described by 
power laws with the same critical exponent. We recover the results of the U(l) slave-rotor formula- 
tion in the presence of a strong breaking of pseudospin symmetry. Our analysis suggests that both 
a third- neighbor electron hopping, and/or pseudospin-breaking terms such as a nearest-neighbor 
density interaction, may help to stabilize a spin liquid phase. 



I. INTRODUCTION 



In the quest for a quantum spin liquid Mott insulator 
in two or more spatial dimensions^ most of the theoret- 
ical attention to date has focused on quantum magnets. 
While this effort has resulted in significant progress in 
our general understanding of a variety of different spin 
liquids, as well as some simple models that have been 
shown to realize such ground states (see Refs. for 
reviews), spin liquids may be rather rare in the simplest 
quantum magnets. In order to advance the understand- 
ing in this area, more models and candidate materials 
are urgently needed, and it is important to look to other 
classes of systems. 

Following a series of recent experiments, the organic 
material k-(BEDT-TTF) 2 Cu 2 (CN) 3 has emerged as the 
first clear candidate for a spin liquid ground state£>£>£ It 
is believed that this material can be described in terms of 
a single-band triangular lattice Hubbard model at half- 
filling (one electron per site). The experimentally ob- 
served spin liquid lies just on the insulating side of a 
metal-insulator transition; this transition and the metal- 
lic state across it, which becomes a superconductor at 
low temperatures, can be accessed by applying pressure. 
The interpretation in terms of the Hubbard model is that 
applying pressure tunes t/U (the ratio of kinetic energy 
to Coulomb repulsion) and drives the system across the 
Mott transition. Among the many striking features ob- 
served thus far, it is notable that the spin liquid physics 



appears adjacent to the Mott transition; this is natu- 
ral, as local charge fluctuations tend to disrupt magnetic 
ordering, which in turn is more likely to be favored for 
stronger correlation (smaller t/U). The message is that, 
when the quantum nature of S = 1/2 spins and geometri- 
cal frustration are not enough to disrupt magnetic order- 
ing, as is likely the case for the triangular lattice Heisen- 
berg model^ local charge fluctuations may be able to tip 
the balance toward a spin liquid state. We are therefore 
optimistic that more examples of spin liquid physics are 
waiting to be found near the Mott transition in other 
systems. 

The primary subject of this paper is the development 
of a general formalism to describe spin liquids adjacent to 
the Mott transition in Hubbard models. At a minimum, 
we wish to describe the simplest scenario, namely a di- 
rect (and continuous) quantum phase transition between 
a spin liquid Mott insulator and a conventional conduct- 
ing state, which may be a metal or a superconductor. 
What features should the necessary formalism possess? 
First, a description starting from the level of mean-field 
theory should be possible. Generally, spin liquids can be 
described in terms of S — 1/2 spinons, which may be 
either bosons or fermions, carrying zero electromagnetic 
charge and coupled to an emergent gauge field. In the 
present case a description in terms of fermionic spinons is 
preferable, as this allows a mean-field description of the 
Mott transition. Therefore, a natural starting point is to 
split the electron into a fermionic spinon and a bosonic 
degree of freedom carrying the charge. 
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For our formalism to be as general as possible, it should 
tie in to the fermionic-spinon theory of the Heisenberg 
model, 1 - which is appropriate for describing a large class 
of spin liquids in the \&vge-U /t limit. In particular, any 
spin liquid state that can be described in the Heisenberg 
model should also be present in our Hubbard model for- 
malism. In all slave-particle formulations a local gauge 
redundancy is introduced; in the fermionic-spinon theory 
of quantum magnets this is an SU(2) redundancy 11 ' 12 In 
order to satisfy the above criterion, our formulation of 
the Hubbard model should also possess an SU(2) gauge 
redundancy. Therefore, we need to generalize the SU(2) 
gauge theory of the Heisenberg model to the Hubbard 
model. This program is in the same spirit as the SU(2) 
gauge theory of the t-J model, which allows for a descrip- 
tion of doped spin liquidsi 13 ' 14 The path pursued here al- 
lows instead for a description of spin liquids with local 
charge fluctuations, as appropriate for insulating phases 
of the half-filled Hubbard model. 

A formalism partially satisfying the above crite- 
ria has already been achieved in the U(l) slave-rotor 
formulation ! 15 i 16 In this approach the electron is formally 
split into a product of a fermionic spinon and a charge- 
carrying rotor boson. There is only a local U(l) redun- 
dancy, however, and only a subset of possible spin liquid 
phases can be described (at the mean-field level) with- 
out extending the formalism; in particular, states with 
spinon pairing and a Zi emergent gauge field, or those 
where the SU(2) gauge symmetry is unbroken, are not 
accessible. 

Here, we build on the U(l) slave-rotor formulation and 
develop the more general SU(2) slave-rotor formulation, 
which satisfies the above criteria for a general descrip- 
tion of spin liquids near the Mott transition. Introduc- 
ing an SU(2) gauge redundancy, the electron is split into 
a product of a fermionic spinon and an "SU(2) matrix 
rotor." While this formalism should be useful quite gen- 
erally, and we do sketch its development for a general 
Hubbard model here, we develop it primarily in the con- 
text of the honeycomb lattice Hubbard model (at half- 
filling), where it turns out to be not only useful but in 
fact essential. This is due to the pseudospin symmetry of 
Hubbard models on bipartite lattices; 17 ' 18 as discussed in 
more detail below, this symmetry is explicitly (and inap- 
propriately) broken by the U(l) slave-rotor formulation, 
but is respected by the SU(2) formulation. This turns 
out to have important consequences for the physics near 
the Mott transition. 

The motivation to consider the honeycomb Hubbard 
model comes largely from the present understanding 
of the k-(BEDT-TTF) 2 Cu 2 (CN) 3 system. Recent the- 
oretical work has proposed that the spin liquid state 
can be described in terms of a Fermi surface of neu- 
tral spin- 1/2 spinons coupled to an emergent U(l) gauge 
fiek h 16 ' 19 ' 20 ' 21 Also, microscopic studies of the triangu- 
lar Hubbard model find a spin- liquid insulator adjacent 
to the Mott transition ) 22 ! 23 although the accuracy of the 
approximations involved is not clear. These ideas are 



quite promising as an explanation of the experimental 
data. However, two obstacles to further progress exist 
in this system, and our understanding would be greatly 
advanced if a complementary system exhibiting similar 
physics but not suffering from these problems can be 
found. First, the spinon Fermi surface state is relatively 
complicated, and our understanding of it is probably no 
so well-developed compared to that of other spin liquid 
states. Second, the triangular lattice Hubbard model 
cannot be effectively studied using quantum Monte Carlo 
simulation due to the notorious sign problem. 

The honeycomb lattice Hubbard model is then an ex- 
cellent candidate to complement the understanding of the 
triangular lattice system. First, because the honeycomb 
lattice is bipartite, the model can be studied by quantum 
Monte Carlo simulation without a sign problem. Also, 
the lack of nesting implies that a metal-insulator transi- 
tion must be present at nonzero U /t, unlike in the square 
lattice. Second, it turns out that the most likely can- 
didate for a spin liquid is an example of an algebraic 
spin liquid (ASL) ; 24 ' 25 ' 26 ' 27 which has only spinon Fermi 
points rather than a full Fermi surface. Such spin liquids 
are a good deal simpler than the spinon Fermi surface 
state, and several recent works have significantly solid- 
ified the theoretical understanding ! 26 ' 27 ' 28 ' 29 ' 30 Further- 
more, it is natural to ask whether a spin liquid is possible 
in the honeycomb model, as it lacks magnetic frustration. 
We speculate that magnetic frustration is not a necessary 
condition for a spin liquid ground state, and that the non- 
trivial competition between kinetic and potential energy 
near the metal-insulator transition can be enough on its 
own. The honeycomb model provides a simple system in 
which to address this issue. At present, we do not know 
of a strongly correlated material described by the hon- 
eycomb Hubbard model; the results of this paper make 
it clear that such materials could be very interesting, if 
they can be fabricated. 

The Hamiltonian of the honeycomb lattice Hubbard 
model is 

w h = - J2 (w4 a <V Q + H.c.) (1) 

(r,r') 

+ U ^(4 T c rT - ^){cl l c rl - -). 

r ~ 

Here, r labels the sites of the honeycomb lattice (shown 
in Fig. [1]), c ra destroys an electron on site r with spin a, 
the sum in the first term is over all pairs of sites (r, r'), 
and repeated spin labels are summed over. We some- 
times find it convenient to label honeycomb lattice sites 
by the pair (R,i), where R is a Bravais lattice vector, 
and i — 1,2 labels the two-site basis of the unit cell 
as shown in Fig. [TJ The honeycomb lattice is bipar- 
tite, and in order to preserve this structure we require 
t rr ' = when r and r' lie in the same sublattice. In 
addition to the nearest-neighbor hopping t, we will also 
discuss the possible effect of the third-neighbor hopping 
t" . We restrict our attention entirely to the case of a 



3 




FIG. 1: The honeycomb lattice with A and B sublattices 
labeled. The unit cell consists of the two sites (labeled 1 
and 2) inside the dashed-line box. Unit cell positions are 
specified by Bravais lattice vectors R = n\a\ + 71202, where 
ni,n2 are integers. The primitive vectors ai = (1,0) and 
a? — (1/2, y/3/2) are shown in the figure. 

half-filled band; combined with the bipartite structure 
this ensures particle-hole symmetry, which is required for 
sign-problem-free Monte Carlo simulation. Our goal here 
is to understand the simplest scenarios for spin liquid 
physics near the Mott transition, and to make predic- 
tions that can be tested in future numerical simulations. 
One quantum Monte Carlo study of the honeycomb Hub- 
bard model has already been carried out^i see Sec. IVIII 
for a discussion. 

This model has a pseudospin SU(2) symmetry, which 
will play a crucial role in our analysis and is reviewed 
in more detail in Sec. IIIIl On a single site of the lat- 
tice, pseudospin rotates between the empty and doubly- 
occupied states, leaving the S = 1/2 singly-occupied 
states unaffected. It will be useful to contemplate adding 
the nearest-neighbor density interaction 

Uv = V (4«<Va - l)(4*aVa " ^ ■ ( 2 ) 
{rr' ) 

This term breaks pseudospin down to an 0(2) subgroup 
generated by U(l) charge rotations and particle-hole 
symmetry. 

The model Eq. (fTJ) has a stable metallic phase for small- 
U/t, and the laxge-U/t ground state is expected to be 
an insulating Neel antiferromagnet (AF). Actually, the 
metallic state is a semimetal (SM) with Dirac nodes in- 
stead of a full Fermi surface, and a linearly- vanishing den- 
sity of states at the Fermi energy. At intermediate-/7/i 
there must therefore be a metal-insulator transition, with 
the possibility of spin liquids or other interesting physics 
nearby. This situation contrasts dramatically with the 
square lattice model, which, at half-filling and in the ab- 
sence of intra-sublattice hopping, is believed to be in the 



(a) 

SM AF 

h X u/t 

SM-AF 

(b) 

SM ASL AF 

| X X " 

SM-ASL ASL-AF 



FIG. 2: Two scenarios for the zero-temperature phase dia- 
gram of the honeycomb Hubbard model as a function of U/t. 
In (a), there is a direct transition between the semimetal (SM) 
and Neel antiferromagnet (AF). This SM-AF transition does 
not involve spin liquid physics. Another possibility, shown in 
(b), is that a region of algebraic spin liquid (ASL) intervenes 
between SM and AF. 



AF phase for all values of U/t, due to the nested Fermi 
surface. 

We caution that the very simplest possibility for the 
metal-insulator transition actually has nothing to do with 
spin liquid physics, and can be understood by coupling 
a fluctuating Neel order-parameter field to the gapless 
Dirac fermions of the semimetal. This results in a direct, 
continuous transition between SM and AF, which can be 
analyzed in a suitable large- N or 4 — e expansion^ The 
phase diagram in this scenario is shown in Fig. [5^. While 
this critical point is already interesting, as the coupling to 
the gapless fermions changes the critical exponents from 
their 0(3) values,— we focus on the even more exciting 
possibility of spin liquid physics near the Mott transition. 

The main result of our analysis is the possibility of a 
Mott transition from the semimetal to an algebraic spin 
liquid Mott insulatori 24 i 25 ' 26 i 27 The resulting phase dia- 
gram is shown in Fig. ^jp. Two transitions are required 
as U/t is increased: first, the Mott transition from the 
semimetal to the ASL, and then a transition from the 
ASL to the Neel antiferromagnet; an example of the lat- 
ter transition has been studied in Ref. [33|. The ASL, 
which is a critical state with power-law correlations for 
many observables and no spin gap, can be thought of 
as a two-dimensional analog of the ground state of the 
S = 1/2 Heisenberg chain. It is an interacting criti- 
cal state in the sense that the fermionic spinons are not 
free quasiparticles even at the lowest energies. It can be 
described by a field theory of massless Dirac fermions 
coupled to a gauge field, and can roughly be thought 
of as the deconfined phase of the gauge theory. Actu- 
ally there are two different kinds of ASL, depending on 
whether the gauge field is U(l) or SU(2) - in the present 
case we find an SU(2) ASL. Both kinds of ASL can be 
analyzed reliably in the limit of a large number of Dirac 
fermion flavors (large- N f limit) and 
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have nontrivial critical exponents that vary as a function 
of Nf. The U(l) ASL is a stable phase in the large- Nf 
limit this is expected also to be true of its SU(2) 
counterpart^ 

For simplicity, in this paper we shall assume that the 
SU(2) ASL is a stable phase. While this is true in the 
large- Nf limit, it need not be the case for the Hubbard 
model. Our results may be relevant to the Hubbard 
model even if the SU(2) ASL is unstable; this is briefly 
discussed in Sec. I VIII 

Our results clarify a brief discussion in Ref. Ilq . which 
focused on the Mott transition in the triangular lattice 
Hubbard model, but also considered the honeycomb lat- 
tice model using the U(l) slave-rotor formulation. That 
analysis missed the role of pseudospin symmetry and 
found a Mott transition to the U(l) ASL. Here, the SU(2) 
slave-rotor formulation has allowed us to properly take 
the pseudospin symmetry into account, and instead we 
find a Mott transition to a SU(2) ASL. This is not just 
a subtle distinction; the U(l) and SU(2) ASL have dis- 
tinct physical properties, and it should be rather easy to 
distinguish them in a numerical simulation. If, however, 
the pseudospin symmetry is broken - by adding Hv , for 
example - a Mott transition to the U(l) ASL again be- 
comes possible. Based on the intuition that gauge the- 
ories with more colors are more susceptible to confine- 
ment, the U(l) ASL is likely to be more stable than its 
SU(2) counterpart. Therefore, in order to observe spin 
liquid physics in numerical simulations it may actually 
be helpful to include Tiy We stress that the issue of 
pseudospin does not affect the results of Ref. for the 
triangular lattice model, where pseudospin symmetry is 
absent. However, it will still be interesting to apply the 
SU(2) formulation there, as it addresses the competition 
among a wider range of spin liquids. 

Another of our major results is a suggestion for how 
to detect the SU(2) ASL in numerical simulations. Very 
recently, Ran and Wen have shown that the SU(2) ASL 
has an emergent Sp(4) symmetry at low energies, which 
contains spin-rotation symmetry as a subgroup^ This 
work was done in the context of the 7r-flux spin liquid on 
the square lattice, which has the same low-energy theory 
as the SU(2) ASL studied here. Furthermore, Ref. Hi 
showed that the Sp(4) symmetry unifies the Neel vec- 
tor and the two-component order parameter for colum- 
nar and box valence-bond solid (VBS) states^ into a 
five-component vector. These observables have slowly- 
decaying power law correlations, all with the same crit- 
ical exponent. These results are analogous to those of 
Ref. [30 on SU(4) emergent symmetry and a unification 
of various "competing orders" in the U(l) ASL. In the 
honeycomb case we also find a unification of the Neel 
vector with the VBS order shown in Fig. [4] Therefore, 
we propose to look for the SU(2) ASL in numerical simu- 
lations by testing for the simultaneous presence of slowly 
decaying Neel and VBS correlations (see Sec. IVIBI for 
further discussion). 

We now give an outline of the paper and sketch our ap- 



proach in some more detail. First, in Sec. HH we briefly 
discuss connections to some earlier work. Following a 
brief discussion of pseudospin symmetry (Sec. lIIIj) . we re- 
view the U(l) slave-rotor formulation in Sec. lIV Al there, 
the electron is split into a fermionic spinon and a charge- 
carrying 0(2) rotor. We show that the U(l) formulation 
explicitly breaks the pseudospin symmetry and is thus 
inadequate to describe the honeycomb Hubbard model. 
We are thus led to introduce the SU(2) slave-rotor formu- 
lation in Sec. lIVBl This is done by splitting the electron 
into fermionic spinons and an SU (2) matrix rotor. In 
contrast to the 0(2) rotor, which can be thought of as 
a quantum particle constrained to move on a ring, this 
rotor is a particle constrained to move in the space of 
SU(2) matrices (or, equivalently, on the 3-sphere). Our 
formulation is obtained by first organizing the electron 
operators into a 2 x 2 matrix 

*, - ( Crl erC h ) . (3) 

where e r = ±1 for r lying in the A/B sublattice 
(see Fig. [T]). This matrix is designed so that left- 
multiplication by an SU(2) matrix is a spin rota- 
tion, while pseudospin rotation is achieved by right- 
multiplication by an SU(2) matrix. We therefore split 
^ r into the product 

*r = F r Z r , (4) 

where F r is a 2 x 2 matrix of fermionic spinon operators 
f ra defined just as in Eq. ([3]), and Z r is an SU(2) matrix. 
Really, Z r is a set of operators giving the position coordi- 
nates of a particle moving in the space of SU(2) matrices. 
There are also corresponding angular momentum opera- 
tors that generate right- and left- rotations of Z r . The 
representation Eq. (|4]) has a local SU(2) redundancy, and, 
upon enforcing an appropriate local constraint, the Hub- 
bard model can be exactly rewritten as an SU(2) gauge 
theory. 

While the development in Sec. IIVBI is specific to the 
Hubbard model on a bipartite lattice, it is straightfor- 
ward to apply this formulation to a general Hubbard 
model with only minor modifications. This is sketched 
in Sec. IIV CI further development of this case is left to 
future work. 

In Sec. [V] a mean-field theory is constructed to 
describe spin liquid and valence-bond solid insulating 
phases. The mean-field theory can also access the Mott 
transition from a given spin liquid phase to a correspond- 
ing conducting ground state. We have considered a re- 
stricted class of mean-field ansatz, in order to focus on 
the simplest possibility of a direct, continuous transition 
from a spin liquid to the semimetal; in Sec. lVBl thc mean- 
field phase diagram is discussed. The saddle point cor- 
responding to the SU(2) ASL is found to be the lowest 
energy state near the Mott transition, and the stability 
of this state is enhanced by the addition of the third- 
neighbor hopping t". 
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In Sec. IVII we develop an effective low-energy theory 
that encapsulates the universal features of the SU(2) 
ASL, as well as those of the Mott transition and the 
conducting state (the semimetal) on the other side. The 
low-energy theory consists of fermionic spinons and rotor 
bosons coupled to a fluctuating, emergent SU(2) gauge 
field. The continuum field theory of the spin liquid con- 
sists of Dirac fermions coupled to the SU(2) gauge field, 
and from this starting point we find enhanced Neel and 
VBS correlations in the spin liquid state. This leads to 
specific suggestions for the detection of the spin liquid in 
numerical simulations (Sec. IVI Bj) . The field theory of the 
Mott transition, which also includes a charge-carrying 
bosonic matter field in the fundamental representation 
of the SU(2) gauge group, is written down and briefly 
discussed in Sec. IVI CI In Sec. IVI D[ the connection be- 
tween pseudospin breaking and the U(l) slave-rotor for- 
mulation is made more concrete. At the level of the effec- 
tive theory, it is shown that pseudospin breaking of the 
type induced by Hv can lead to a phase transition from 
the SU(2) ASL to the U(l) ASL. This spin liquid state, 
as well as its own direct transition to the semimetal, are 
described by the effective U(l) gauge theory of Ref. [H, 
which therefore applies here in the presence of strong 
enough pseudospin breaking. 

In Sec. lVIll we summarize our main results for the hon- 
eycomb lattice model, with a focus on implications for 
future numerical simulations. Furthermore, we briefly 
discuss the possibilities for applying the SU(2) slave- 
rotor formulation to more general Hubbard models, and 
also speculate on its possible utility in obtaining a better 
physical picture of fermionic spinons. Various technical 
details are described in the appendices. 



II. RELATION TO EARLIER WORK 

The SU(2) slave-rotor formulation has recently been 
constructed independently byK.-S. Kim via a functional 
integral approach. 40 ' 41 Refs. iolfill began with the square 
lattice Hubbard model and decoupled the on-site repul- 
sive interaction with a pseudospin vector field. Then, af- 
ter a series of changes of variables, the microscopic SU(2) 
slave-rotor functional integral of Eqs. (|45ti47|) was ob- 
tained. Here, in contrast, we introduce the SU(2) slave- 
rotor variables working in the Hamiltonian, and then go 
to a functional integral. Our approach makes clear that 
this formulation is a generalization of the SU(2) gauge 
theory of the Heisenberg model to the case of the Hub- 
bard model. 

Despite the equivalence of our formulation to that of 
Kim, we believe that this formulation has been ana- 
lyzed incorrectly in Refs. l40H4ll . Upon moving away from 
half-filling, Refs. I40ll4ll found a non-Fermi liquid metal- 
lic phase intervening between a conventional Fermi liq- 
uid and a spin liquid state. First, the spin liquid state 
certainly cannot persist away from half-filling; the ro- 
tor bosons will enter the ground state at finite density 



and condense, resulting in a conducting ground state (ei- 
ther a Fermi liquid metal or a superconductor) with the 
gauge field in the Higgs phase. Second, the SU(2) slave- 
rotor formulation probably cannot be used to describe 
a non-Fermi liquid metallic phase. We believe the latter 
error arises because Refs. l4Gl4ll do not correctly take into 
account the separate roles of pseudospin symmetry and 
SU(2) gauge redundancy. We comment on this in more 
detail at the end of Sec. IVI Al 

We also mention earlier work by Wu et. at, where the 
idea of quantum disordering a <i-wave superconductor by 
working with the 0(3)-vector pseudospin order param- 
eter was explored^ These ideas are connected to the 
approach discussed here. 



III. REVIEW OF PSEUDOSPIN 

Here we briefly review the pseudospin symmetry of 
Hubbard models on bipartite lattices. 17 ' 18 While pseu- 
dospin symmetry will not be present (except perhaps ap- 
proximately) in any realistic Hubbard model, even if the 
lattice is bipartite, it is still important to consider it for 
the analysis in this paper. First, it is present for the 
simplest Hubbard models on bipartite lattices, namely 
for Hamiltonians with only bipartite hopping and on-site 
repulsion. Because we propose to study the honeycomb 
Hubbard model via quantum Monte Carlo simulation, it 
is important to take pseudospin into account. Second, it 
will turn out to lead us to a generalization of the SU(2) 
gauge theory of the Heisenberg model to the case of the 
Hubbard model, which will be useful much more gener- 
ally 

The most familiar symmetries of the Hubbard Hamil- 
tonian Eq. |T|) , in addition to particle-hole symmetry, are 
the space group symmetry, time reversal, U(l) charge 
rotations generated by the total electron number, and 
SU(2) spin rotation generated by Stot = J2 r where 
Sr = c la (T af3C r p/2 and er = (a 1 , a 2 , a 3 ) is a vector of the 
2x2 Pauli matrices. The pseudospin symmetry is most 
easily exposed by organizing the electron operators into 
the 2x2 matrix ^ r as in Eq. Left-SU(2) rotations 
of ^ r are generated by the spin density: 



„iaS, 



-ict-S, 



(5) 



The pseudospin density T r is a spin singlet and is de- 
fined as 



T z = c ) 

T r € r C r ^C r ^ 



where T± = T* ± iT$, and 



1, r e A 
-1, r e B 



(G) 

(7) 
(8) 



(9) 
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The pseudospin density generates right-SU(2) rotations 
of * r : 

e laTr y r e- iaT T = ty r e iarT/2 . (10) 

The Hubbard Hamiltonian is manifestly invariant both 
under spin and pseudospin symmetry in the form 

«H = - E ^tr($tf r ,) + y^(T r 2 ^), (11) 

(r,r') r 

where the first term is already Hermitian without taking 
the complex conjugate. It should be noted that the bi- 
partite nature of the lattice was crucial in obtaining this 
form; it entered via the factors of e r in the definitions 
of $ r and T r . It is also useful to note that spin and 
pseudospin densities can be expressed simply in terms of 

S r = |tr(*J,o-* r ) (12) 

Tr = |tr(\£r r <7*t). (13) 

Pseudospin is an enlargement of the 0(2) symmetry 
group generated by U(l) charge rotations and particle- 
hole symmetry. Under a U(l) charge rotation c ra — > 

e~ l9 c ra and hence $ r — > ^ r e~ %ea , so the conserved 
U(l) charge density is simply — 2T*, as is apparent from 
Eq. ([6]) . Particle-hole symmetry acts on the electron op- 
erators as c ra — > e r (icr 2 ) a pcl,p, and so ty r — > ^f r (— i<J 2 )- 
Particle-hole symmetry is therefore implemented by the 
unitary operator exp (— Ylr ^r) > which is again part 
of the pseudospin symmetry group. 

To understand a little better the physical meaning of 
pseudospin symmetry, it is useful to consider the stag- 
gered pseudospin = e r T r . The z-component (T s ) z 
is the order parameter for a staggered charge density 
wave state, while the transverse components make up 
the order parameter for an s-wave superconductor. The 
pseudospin symmetry rotates among these superficially 
unrelated states. 

Pseudospin is broken by the nearest-neighbor repulsion 
term Tiy of Eq. @; in terms of pseudospin generators 
TLy can be written 

Hy=4l/^T r 2 T;„ (14) 

(rr') 

This is an Ising exchange in pseudospin space. It breaks 
the pseudospin symmetry down to the 0(2) subgroup 
generated by U(l) charge rotations and particle-hole 
symmetry. 

IV. SLAVE-PARTICLE FORMULATION 

Here, we develop the SU(2) slave-rotor formulation. 
We are led to this formulation in the context of Hubbard 



models on bipartite lattices (e.g. the honeycomb lattice), 
as a means to properly account for the pseudospin sym- 
metry. However, its real power is as an extension of the 
SU(2) gauge theory of the Heisenberg model to the case 
of the Hubbard model. As such, it should be useful much 
more generally, even in the absence of pseudospin sym- 
metry. 

First, in Sec. II V A[ we review the U(l) slave-rotor 
formulation ; 15 i 16 and explain why it is not adequate to 
describe the honeycomb lattice model. Next we build on 
the U(l) formulation and develop the SU(2) slave-rotor 
formulation in the context of a Hubbard model on a bi- 
partite lattice (Sec lIVB]) . Finally, in Scc. lIVCl we sketch 
the development of the SU(2) formulation for a general 
Hubbard model, without the assumption of a bipartite 
lattice. The remainder of the paper is devoted to an anal- 
ysis of the bipartite case and specifically the honeycomb 
lattice model. Applications of the general formalism are 
left to future work. 



A. Attempt at a U(l) formulation 

The U(l) formulation begins by rewriting the electron 
as the product of a charge-neutral fermionic spinon and 
a charge-carrying 0(2) rotor 

This is a rewriting of the electron in terms of spin-charge 
separated variables. In Eq. (|15|) . f ra is the spinon de- 
struction operator, and the operator e l ^ r measures the 
position coordinate of the 0(2) rotor, which is con- 
strained to lie on the unit circle in the complex plane. 
There is also a conjugate rotor angular momentum L r 
taking integer eigenvalues, and satisfying the commu- 
tation relation [L r ,e ±l< ^'] = ±5 rr , e ±l ^ . The Hilbert 
space on each site is a product of spinon and rotor Hilbert 
spaces, and a local constraint must be imposed to restrict 
this to the physical Hilbert space of the original Hubbard 
model. This is accomplished by requiring 

fUra -L r = l. (16) 

In these variables there is a local U(l) redundancy, where 
all physical operators are invariant under the transforma- 
tion f ra — > e iXr f ra and 4> r — > <f> r + X r . 

As it stands, Eqs. fT5|) and (fl~6|) can be used to exactly 
rewrite the Hubbard Hamiltonian as a U(l) gauge the- 
ory, but this does not yet provide any information about 
the low-energy physics. To proceed, one can develop a 
mean-field theory in terms of the fermions and rotors and 
expand in the fluctuations about the mean field saddle 
point. This was done for the honeycomb Hubbard model 
in Ref. [lri the result can be encapsulated in a Hamil- 
tonian U(l) lattice gauge theory, and in the mean-field 
analysis of the honeycomb Hubbard model a U(l) alge- 
braic spin liquid was found adjacent to the Mott transi- 
tion. In the low-energy Hamiltonian, the local constraint 
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is modified to 

(dive) r = f r J ra -L r -1, (17) 

where e rr > is the integer- valued emergent electric field ly- 
ing on the links of the honeycomb lattice, and (div e) r — 
J2 r > n n r e rr' is its lattice divergence. There is a term 
in the effective Hamiltonian of the form h ^2/ rr /\ e 2 r , > so : 
when h < oo, (e 2 r ,) ^ and the local constraint has 
been softened. For ft, = oo and a hard constraint, the 
only states of the rotor angular momentum allowed are 
L = 0, ±1. However, in order to construct a useful effec- 
tive theory the constraint must be softened, and higher 
angular momentum states will be present. 

In terms of pseudospin, the L = state is a pseudospin 
singlet, and the L = ±1 states form a pseudospin doublet. 
But the higher angular momentum states do not trans- 
form as irreducible representations of the pseudospin, and 
as soon as they are present, the resulting effective Hamil- 
tonian explicitly breaks the pseudospin symmetry. The 
U(l) formulation therefore leads to effective theories that 
do not respect the pseudospin symmetry; this invalidates 
its results, at least near the metal-insulator transition. 

If the pseudospin symmetry is broken explicitly in the 
microscopic Hamiltonian - for example, by adding the 
term Hy of Eq. (J2J) - the U(l) formulation remains valid, 
in the sense that it leads to a sensible low-energy theory 
respecting the symmetries of the model. The low-energy 
theories obtained in the U(l) formulation can be invari- 
ant under the 0(2) subgroup of the pseudospin generated 
by U(l) charge rotations and particle-hole symmetry, so 
the formulation can be used in models with these symme- 
tries, or with less symmetry (as in the triangular lattice 
Hubbard model, which was the focus of Ref. Ha)^ 

In fact, we can work backwards from the effective U(l) 
lattice gauge theory Hamiltonian to obtain a modified 
Hubbard model, which itself breaks pseudospin symme- 
try. This is done by taking h large compared to all other 
energies in the effective Hamiltonian, and doing degener- 
ate perturbation theory in the manifold of states where 
e rr > = 0. One term that will be present is a spinon hop- 
ping term with coefficient t s . At order t s /h in perturba- 
tion theory (the leading order), two terms are generated. 
The first is an antiferromagnetic Heisenberg exchange. 
The second is nothing but the repulsive interaction Hv, 
with V oc tg/h. The interaction so obtained is repulsive, 
and it is natural to guess that positive (rather than neg- 
ative) V is more likely to favor a U(l) spin liquid state. 

B. SU(2) slave-rotor formulation and 
spin-pseudospin separation 

The discussion above shows that the U(l) formula- 
tion needs to be generalized to respect the pseudospin 
symmetry. To accomplish this, we need to generalize 
spin-charge separation to spin-pseudospin separation - 
recall that the charge-rotation symmetry is a subgroup 
of the pseudospin. As is apparent from the ^ r matrix of 



Eq. Q , the electron operator transforms as a doublet un- 
der both spin and pseudospin. Therefore we should split 
it into a doublet of fcrmionic spinons, and, in contrast to 
the U(l) formulation, a doublet of bosonic objects carry- 
ing the pseudospin. As in the U(l) formulation, we shall 
take this bosonic object to be a kind of rotor. 

We proceed in this manner to develop the SU(2) slave- 
rotor formulation. We begin by writing 

m r = F r Z r , (18) 

where 

Fr = ( ^ e _ rf \\ ) (19) 
V Jrl £rj r ^ J 

is a matrix of S = 1/2 spinon creation and destruction 
operators with anticommutation relations {fra, flip} — 
8 a p8 rr > , and Z r is an SU(2) matrix. Actually Z r is a 2 x 2 
matrix of operators measuring the position coordinate 
of an SU(2)-matrix rotor, which can be thought of as 
a quantum particle constrained to move in the space of 
SU(2) matrices. The matrix Z r can be parametrized by 
writing 

Zr = ( z :i J 2 ), (20) 

\ z r2 z rl J 

where the z r , are operators, z* ri = z' ri , and there is a con- 
straint |z r i| 2 + \z r 2\ 2 = 1. This parametrization makes 
clear that the SU(2) matrix rotor is identical to the more 
familiar 0(4) rotor, a quantum particle constrained to 
the 3-sphere. Furthermore, it is important to note that 
the decomposition Eq. (|18p contains the minimal num- 
ber of degrees of freedom for spin-pseudospin separa- 
tion: a single doublet of spinons, and a single doublet 
of pseudospin-carrying rotor bosons. 

There is a local SU(2) redundancy, and a corresponding 
invariance under the transformation F r F r e~ laT ' a I 2 
and Z r — > e lOLr ' rT /' 2 Z r . This will allow us to rewrite the 
Hubbard model exactly as an SU(2) gauge theory. In 
fact, focusing only on the spinons, as is appropriate in 
the U — > oo limit, the SU(2) gauge structure reduces to 
the one present in the Heisenberg model . 11 ' 12 This was 
extended to to the case of a doped system described by 
the t-J model in Refs. [l3| and In the present work, 
we show that this SU(2) gauge structure extends to the 
full Hubbard model. 

The enlarged slave-particle Hilbert space on each site 
consists of a product of the Hilbert space of S = 1/2 
fcrmionic spinons and that of the SU(2) matrix rotor. To 
understand this in more detail we first state some facts 
about the SU(2) matrix rotor that are explained more 
completely in Appendix[A] The Hilbert space for a single 
SU(2) matrix rotor has a basis of position eigenstates 
{\Z)} satisfying Z\Z) = Z\Z)\ here, Z e SU(2) and Z 
is the 2x2 matrix of operators discussed above, and 
we use the "hat" symbol as a reminder that an object 
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is an operator when it is not otherwise clear from the 
context. We define the analog of angular momentum 
operators Jr and Jl as the generators of right- and left- 
SU(2) rotations: 



e ia - jR \Z) 
e ta - jL \Z) 



iZe-™-"/ 2 ), 
|e la - CT/2 Z). 

This implies the operator identities 

ux-Jl gg— tot- Jl _ e —ia-tr/2g 

and the commutation relations 



(21) 
(22) 



(23) 
(24) 





= ie ijk 4 


(25) 


[4,4] 


= ie ijk J% 


(26) 




= 


(27) 




-*(t) 


(28) 


[4,z] 


- -(t)*- 


(29) 



We can also work with the basis of angular momen- 
tum eigenstates specified by {\£r =£,£l = £,rnR,mL}}, 
where £r and II are the total right and left angular 
momenta, respectively, and thr and are the corre- 
sponding projections along the z-axis. It is shown in 
Appendix [5] that only states with Ir — II are allowed 
in the Hilbert space; this is a higher-dimensional ana- 
log of the more familiar fact that, for a particle moving 
on the 2-sphere, only integer angular momenta are al- 
lowed. The unique rotationally invariant state is denoted 
|0) rot = |0,0,0,0>. 

Now we can return to the Hubbard model, where we 
place an SU(2) matrix rotor on every site. The angular 
momenta are denoted by Jr{v) and Jz,(r). The pseu- 
dospin is clearly identified with Jr(t): 



Jr(t). 



(30) 



We now have enough information to represent the phys- 
ical states for one site r of the Hubbard in the enlarged 
slave-particle Hilbert space. The physical vacuum state 
for the electrons (the empty state) is denoted by |0), and, 
in the slave-particle Hilbert space, we define a vacuum 
state |0) sp = |0)/ <g> |0) ro t, where |0)/ is the spinon vac- 
uum. Then we have 



4 T I°) «-> /J t i°> s p 
4 jo) 

10) 

4 T 4jo) 



(31) 

/. T J0) SP (32) 
(z rl +e r z; 2 / r t T / r t i )|0) sp (33) 
{<ifUftl - e^2)|0) sp . (34) 



It can be shown that the electron operator in the slave- 
particle Hilbert space, as defined by Eq. (fT5)) . has the 



correct matrix elements among all these states, as well 
as the correct anti-commutation relations. 

We still need a local constraint in order to eliminate 
the unphysical states in the slave-particle Hilbert space. 
As in other slave-particle gauge theories, the appropriate 
constraint is related to the SU(2) redundancy; in fact, it 
involves the generator of SU(2) gauge transformations. 
We define 

Jjr(r) = itr(F P< rFt) > (35) 

which generates right-SU(2) rotations of F r , and 

J G (r) = J F (r) + J L (r). (36) 

We have 

e i a -Ja(r) Zre -ia-J G (r) = g-ia-o-/2^ ^ 

e ia-J G (r) p^ e -ia-J G (r) _ ^ j-gg-j 

so that Ja( r ) is the generator of SU(2) gauge transfor- 
mations; equivalently, it is the SU(2) gauge charge. It 
can be shown that, under the local constraint 



Ja(r) = 0, 



(39) 



the slave-particle Hilbert space is identical to the physical 
Hilbert space of the Hubbard model. 

The relationship of this formulation to the SU(2) gauge 
theory of the Heiscnberg model is more transparent if we 
introduce the spinor 



tr 



(40) 



which transforms as a doublet under SU(2) gauge rota- 
tion. In terms of i/v the fermion part of the SU(2) gauge 
charge takes the more familiar form 



(41) 



Now we can construct the functional integral repre- 
sentation of the Hubbard model in the SU(2) slave rotor 
formulation, using the results above, and those of Ap- 
pendix |XJ Dropping overall additive constants, we write 
the Hamiltonian as 



(r,r') 



with the constraint 



J F {r)+ J L (r) = 0. 



(42) 
(43) 



We can follow the methods of Appendix [Alto write down 
the functional integral 



Z = / VZVFVa T e 



(44) 
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where S = Sz + Sf + S t . We concentrate on zero- 
temperature; the imaginary time r ranges from — oo to oo 
and we write f dr = f°° dr. We have 



(45) 

where the symbol d T is a derivative acting on the left. 
Also, 

Sf = y~) dr f ra d T f ra +ia T -J F (r,T) 



/ tr F,.(5 T 



0*2 



and 



J dr Y t rr ,tr ({F r Z r y (F r > Z r >)) . 



(46) 



(47) 



(r ,r' ) 



So far, this is an exact rewriting of the Hubbard model 
as a (strong coupling) SU(2) gauge theory. Both spin and 
pseudospin symmetry arc manifest in this representation. 
A general gauge transformation acts as follows: 

Z r ( T ) -» Gi(T)Z r (r) (48) 
Fr(r) -> F r {T)G r {r) (49) 

^ - Gl(r)(^ + T )G r (T), (50) 
where G r (r) is an SU(2) matrix. 



C. SU(2) slave-rotor formulation (general lattice) 

The SU(2) slave-rotor formulation is really a gener- 
alization of the SU(2) gauge theory of the Heisenberg 
model to any Hubbard model. In particular, pseudospin 
symmetry is not required, even though we have been fo- 
cussing on models that possess it up to now. Here, we 
digress from the honeycomb lattice problem and develop 
the SU(2) slave-rotor formalism for a general Hubbard 
model, without the assumption of a bipartite lattice and 
pseudospin symmetry; we also include an arbitrary chem- 
ical potential. It turns out that only very simple mod- 
ifications are required. Because we shall not return to 
this more general case later in the paper, it is most eco- 
nomical to proceed by changing notation. This section 
is therefore meant to be self-contained, and the notation 
here is not consistent with the rest of the paper. 

We organize the electron operators into the the 2x2 
matrix 



(51) 



•>1 



The difference from Eq. ([3]) is the absence of e r , which 
is not defined for a general lattice. Spin rotations act on 
^ r as before 



and, also as before, U(l) charge rotations act by sending 
^ r — > ^> r e~ l8a . Even though pseudospin is not in gen- 
eral a symmetry, it is useful to note that the generators 
of right-SU(2) rotations of ^> r are T r , so that, 



where 



and 



T z = -fl-c t c ) 

v 2 ra^ra J i 



(53) 



(54) 



(55) 



Also, T~ = (T+)t and T± = T* ± %TJ>. Dropping addi- 
tive constants and including a chemical potential /Lt, the 
Hubbard Hamiltonian can now be written 



(r,r') r 

(56) 

where the difference from Eq. (fTTj) is the presence of the 
cr matrix in the hopping term. 
We make the decomposition 



where 



ty r — F r Z r , 



rp _ ( frf fL \ 
V frl -frj' 



(57) 



(58) 



and Z r is defined as before, ty r is invariant under the 
SU(2) gauge transformation F r — ► F r e ta ' r ' cr and Z r — ► 
g-ja T .-<7- Again, we define the generators of right and 
left rotations of Z r to be Jr{t) and Jr,(r), respectively, 
and we identify Jr(t) = T r . Defining 



J F (r) = itr (F r trF r t), 

the SU(2) gauge charge is 

J G (r) = J F {r)+J L {r) 



(59) 



(60) 



The local constraint Jg{t) = restricts the slave-particle 
Hilbert space to the physical Hilbert space of the Hub- 
bard model. 

Following the methods of Appendix [A] we can write 
down the functional integral. The partition function is 



Z = / VZVFVa T e- 



(61) 



where S = Sz + Sf + St, Sf is identical to Eq. (|4"B"|) upon 
replacing F r there by the definition Eq. and 
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Sz = W- 



drtr 



zt{d T - 



ia T 



)Z r 



3/i 
2U 



— > drtr 



ia r 



)Z r 



(62) 



Finally, we have 



J dr f rr> tr (a z {F r Z r ^{F r ,Z rl )). 



(63) 



ir,r') 



We have arrived at a representation of a general Hub- 
bard model as an SU(2) gauge theory. The rest of this 
paper focuses on the formulation specific to bipartite lat- 
tices developed in Sec. IIVBI Using the following anal- 
ysis as a guide, it should be possible in future work to 
use the more general formulation to investigate spin liq- 
uid physics near the Mott transition in a panoply of 
models. Motivated by the experiments on k-(BEDT- 
TTF)2Cu2(CN)3, it would be particularly interesting to 
consider the triangular lattice Hubbard model and com- 
pare to the results of Ref. [l6|. 



V. MEAN-FIELD THEORY 

We formulate a mean-field theory based on the gauge 
theory action for bipartite Hubbard models described 
above. This will allow us to access spin-liquid Mott insu- 
lators, valence-bond solid states, and conducting states 
where the charge-carrying rotor bosons have condensed. 
First, we develop the general structure of the mean- field 
theory (Sec. IV Af . Next, focusing on the honeycomb lat- 
tice, we consider a restricted class of mean-field ansatz 
and discuss the resulting phase diagram (Sec. IV Bj) . 



A. General structure of the mean-field theory 

The starting point for the mean-field theory is the func- 
tional integral derived in Sec. IIV Bl We begin by trading 
in the constraint Z r {r) e SU(2) for a Lagrange multi- 
plier. The simplest way to do this is to parametrize Z in 
terms of z\ and zi as in Eq. (|20[) . allowing z\ and Z2 to be 
arbitrary complex numbers, and then imposing the con- 
straint \z\\ 2 + \ z2\ 2 = 1 with a single Lagrange multiplier 
\ r {j). This results in a contribution to the action 



S\ = i dr A r (r) 



itr(2jZ„)-l 



(64) 



To decouple the hopping term St, we follow Ref. [if 
and use the identity 



- / driijdrj-jexp -e(\r] l3 ;| -r^-o^ -r)!jPij) 



(65) 



We write 



St 



(r,T* ; ) 



,{Z r ,Zl) a p{FlF r ,) Pc 



(66) 



Putting e = |i rl .'|AT and performing the decoupling us- 
ing the complex fields ifjp , we trade St for the following 
contributions to the action: 



Sri — 



Stz = 
StF = 



dr iwMorv* 



(r,r ; ) 



JdrJ2 \Ur>\tr {Zt V rr ' Z rl ) 



(r,T*') 



dr t rr ,tv(F r ,(r, rr ')iF}). 

(r,r') 



(67) 
(68) 
(69) 



By decoupling in this manner, we give each bond an ori- 
entation. Because we are working with a bipartite lattice, 
it is natural to take r in the A sublattice and r' in the B 
sublattice in the equations above. This orientation disap- 
pears in many of the mean-field saddle points discussed 
below. 

We shall now parametrize r] rr in a more conve- 
nient form. Because non-invertible matrices make up a 
measure-zero subset of all matrices, we may take T] rr to 
be invertiblc. Then rf r can be expressed as a complex 
number times an arbitrary matrix with unit determinant 
(such matrices form the group SL 2 (C)). In fact, up to a 
set of measure zero, if M £ SL2(C) then M = exp(ibi(j l ) 
for some set of complex numbers 6^.— Therefore, we can 
write 



rr 

V = 



where \xrr' 

\Xrr'\ > 0. 



Xr 



exp 



(70) 



cT r and d^ r are real numbers, and 



Eventually we will look for saddle points where we in- 
tegrate out Z r and F r and work with a free energy de- 
pending on A, a T and r) rr . We want to understand what 
kind of rj r " r matrices give Hermitian Stz and StF, and 
thus a real free energy. Making use of Eq. (|4"0"|) we have 

tr (F^ (v rr 'yF r ) = 4>l {rf r ')H r + ^ r a 2 (v rr ')*^ r ■ 

(71) 

Using Eq. ([70]) . we see that StF is Hermitian if and only 
if 8 rr r = and d|" r = 0, so that rf r is a real number 
times an SU(2) matrix. To proceed similarly with Stz, 
we define 



Z r i 



(72) 
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and we have 



tr (Zl V rr 'Z r ,) = Cv rr 'C r/ + Cl* 2 (v rr ') T <T 2 ( r . (73) 

Again, this is Hermitian if and only if rf r is a real num- 
ber times an SU(2) matrix. 

We now write down the saddle-point equations: 



1 = l -{tr{ZlZ r )) 

= ^ja T + ^j1m(tr(d T Zt*Z r )) 



(v rr ') f = (z r ,zi) 

\t rr l I \ 



(74) 

(75) 
(76) 
(77) 



Following Ref . [l6|, it is useful to make an analytic con- 
tinuation in order to access a broader class of saddle 
points. We consider 



and therefore 

_ / 

r r 

V 



d rr ' 

a T (r) 



\Xrr-W. 



t\ r 

idf 

id T (r), 



(78) 
(79) 
(80) 



rr 

Z 



(82) 
(83) 



where \x z rr ,\ = {XrAe^ , = \Xrr'\e^', U r / 

exp[i(c rr ' ~ d?)a l ], and U r F r = {U F r 'Y = exp [ - 

i(c*" r + d\ r )cr 1 ]. The saddle-point free energy is now 
real, and the saddle-point equations become 



\Xrr>W. 
\X Z rAU T Z 



r r 



- (tr(ZtZ r )) 

^jar^ + ^jlm^xidrZtaZ,)) 

- {ll>t(Tll) r ) 

z„,zl 



F l F r' 

|i rr / 1 



(84) 

> 

(85) 
(86) 
(87) 



These expectation values, and also the mean-field 
ground state energy (see Appendix [Cj , are evaluated us- 



ing the saddle-point action S — Sq + Sf + Sz, where 

So = j dr{ £ M\xU\x z r Muyu r /) 

(r.r') 

r 

Sf = f rfr | E ^<9rVv -^E^ [M*0 ■ er] tp r 
J r r 

+ E W|xfr'l(^^F r >w+H.c.)} (89) 

{z2£[ijl- d2 r+(«r(r)-<r)dr]+K 



(r,r' ) 

S z = I dr 



£ \U r '\\x Z r'\(CU r z r 'Cr'+C.C.)}. 



Cr- 
(90) 



(r ,r' ) 



Mean-field phase diagram and relation to U(l) 
calculation 



Now we focus specifically on the honeycomb Hubbard 
model at half filling, where we are interested in describing 
spin-liquid states adjacent to the Mott transition to the 
semimetal. We restrict our attention to an ansatz of the 
form 



U r z r = 



U r F r 



e- ra * 

e *C' 
e~ r < 



a T {r) = 0. 



(91) 

(92) 
(93) 



I Xrr' 



and | Xrr 'I are arbitrary. While it might be inter- 
esting to determine the mean-field phase diagram consid- 
ering a completely general ansatz, we have not done this. 
First, we are interested in the possibility of a continuous 
transition from the semimetal to a spin liquid Mott insu- 
lator, and as we shall see the above ansatz allows us to 
describe the simplest scenario in which this can happen. 
Second, the mean-field theory only gives very crude infor- 
mation about the energetics of various competing states, 
and is probably only useful as a rough guide to numeri- 
cal simulations. Therefore the value of a more complete 
mean-field study is rather dubious, at least until more 
numerical data are available. 

For the ansatz considered, it turns out that the mean- 
field calculation is equivalent to that already carried out 
for the U(l) slave rotor formulation. 16 This equivalence is 
demonstrated in Appendix [Cl the discussion there shows 
that the results from the U(l) calculation can be used 
here with no change except a rescaling of the onsite inter- 
action by U — > (3/4)£7. It should be stressed that, while 
the mean-field calculations in the U(l) and SU(2) formu- 
lations are formally equivalent (for this restricted class 
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of ansatz), the interpretation of the resulting phase dia- 
gram in the two cases is different. In particular, the spin 
liquid state we find here is described by Dirac fermions 
coupled to an SU(2) gauge field, not a U(l) gauge field 
as was found in Ref. Ilq . 

Now we discuss the mean-field phase diagram, us- 
ing the results of Ref. We consider both nearest- 
neighbor (t) and third-neighbor (t") hopping. For t" = 0, 
the SU(2) ASL is the lowest-energy mean-field state for 
1.26 < U/t < 1.30. For U/t < 1.26 the semimetal is the 
ground state; this can be determined by noticing that 
the gap to Z-boson excitations vanishes as U/t is taken 
to 1.26 from above. For U/t > 1.30, a valence-bond solid 
is the ground state. The third neighbor hopping helps 
stabilize the spin liquid; for t" = — 0.4t, the window of 
SU(2) ASL expands to 1.7 < U/t < 1.9. 

In the case t" = 0, the mean-field critical value 
(U/t) c = 1.26 is significantly less than the value (U/t) c ~ 
4—5 determined by quantum Monte Carlo simulation and 
series expansion^ Of course, it is not expected that the 
mean field theory will be quantitatively accurate, and it 
may not be qualitatively accurate. However, it is rea- 
sonable to use the mean-field results as a rough guide 
to the energetics of the spin liquid state and some of 
its competitors, so it is a useful conclusion that adding 
third-neighbor hopping is likely to enhance the stability 
of the spin liquid. 

We note that for \t"/t\ large enough, the band struc- 
ture changes qualitatively from the t" = case. This can 
be understood by considering the case t — 0, t" ^ 0, 
where the honeycomb lattice breaks into four larger hon- 
eycomb lattices, leading to eight Dirac nodes at the Fermi 
level as opposed to two when t 0, t" ^ 0. Straight- 
forward calculation shows that the two-node band struc- 
ture of t" = is confined to the region —1 < t" /t < 1/3, 
with eight nodes present elsewhere. The free-electron 
Hamiltonian with t" ft = —0.4 lies safely in the two-node 
region, but this is not a priori obvious for the mean- 
field spin liquid state obtained for the same parameters. 
However, it turns out that, for these parameters, the ra- 
tio of spinon hopping matrix elements lies in the range 
—0.4 < t'//t s < 0^ and the spin liquid state thus has 
only two nodes. 

Although we have not pursued this here, it would be 
interesting to investigate the possibility of spin liquid 
physics at larger \t" /t\. One might be able to obtain 
a SU(2) ASL with Nf — 8 Dirac fermions, which should 
be substantially more stable than the Nf = 2 case found 
here (at least as far as local stability to weak perturba- 
tions). The semimetal state will also be different in this 
regime, with eight nodes instead of two. 



VI. LOW-ENERGY EFFECTIVE THEORY 

We shall write down a low-energy theory for the hon- 
eycomb Hubbard model describing the SU(2) ASL, the 
semimetal and the Mott transition between them, as well 



as other nearby phases and critical points. This effective 
theory contains the fermionic spinons and rotor bosons 
coupled to an SU(2) gauge field. Working on a space- 
time lattice, the action is written down in Sec. IVI Al and 
the phase diagram is discussed. Next, in Sec. IVI Bl the 
continuum quantum field theory for the ASL phase is ob- 
tained, and some of its physical properties are discussed. 
The continuum theory of the Mott transition itself is con- 
structed in Sec. IVI Cl Finally, in Sec. IVI Dl we consider 
the effect of breaking the pseudospin symmetry down 
to the 0(2) subgroup generated by particle- hole sym- 
metry and U(l) charge rotations. For sufficiently strong 
pseudospin-breaking of this type, we find a phase transi- 
tion from the SU(2) ASL to a U(l) ASL, and recover the 
results of the U(l) slave-rotor formulation. 

A. Gauge theory description of spin liquid and 
semimetal 

Here we construct the low-energy effective theory and 
use it to give a description of the spin liquid and 
semimetal phases. In principle we could do this, starting 
in the spin liquid phase, by doing a detailed expansion 
about the mean-field saddle point. However, because we 
expect that all fluctuations about this saddle point will 
be massive, except for those captured as an SU(2) gauge 
boson, we can simply write down an SU(2) gauge theory. 
It is convenient to discretize imaginary time and work on 
a space-time lattice, going to a continuum formulation 
later on. 

We consider a Euclidean space-time of stacked honey- 
comb lattices, with lattice sites labeled by by x = (r, r) = 
r + tz, where r = en is a discrete imaginary time (n is 
an integer and e is the time step), r labels the sites of the 
honeycomb lattice, and z is a unit vector pointing along 
the imaginary time direction. On each lattice point we 
have an SU(2) matrix Z x and the 2x2 matrix of fermions 

F x = ( f ;f\ e 4 i( ;\ ) ■ (94) 

On nearest-neighbor bonds within each honeycomb layer 
("spacelike" bonds), as well as on all "timclike" bonds 
connecting two adjacent layers, we place the SU(2) gauge 
field U xx i. The action takes the form S cS = Sf + Sf + 
Sg S , and the partition function is 

Z cB = / VZVfVfVUe' s ° t! . (95) 

Here, 

= C tr {Fx+exUx+ez,xFx) 

X 

+ cEE tr ( f ^^)' ( 96 ) 

t (rr f ) 

where it should be noted that the first term defines an 
orientation on the timelike bonds; that is, the overall 
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sign is reversed if x and x + ez are interchanged. It 
is worth noting that this does not break time-reversal 
symmetry, as time-reversal maps the imaginary time r 
to itself, not to — r. The second term, consistent with the 
spatial reflection symmetries of the honeycomb lattice, 
does not introduce an orientation for the spacelikc bonds. 
Next we have 



ccff 



"*c tr (zl +tz U x+ , z 



t (rr') 



and 



ceff _ 
a g - 



-K 



n 



tr (U XX >U X 



.) +c.c), 



(97) 



(98) 



where the sum is over two types of lattice plaquettes: 
hexagonal plaquettes lying within the honeycomb layers, 
and 4-sided plaquettes containing two timelike and two 
spacelike bonds. The trace is over a product of U xx >, 
which is taken moving along the plaquette edge in an 
arbitrary direction. 

Upon taking the time-continuum limit e — > 0, we re- 
cover the correct dynamics for the fermions and rotor 
bosons. We focus on the terms involving timelike bonds 
Letting t T s = -1/2, t T c = 3/(2Ue) and 



in Sf and Sf 



( ia T - a 
U x+ez<x = exp - e — - — 



(99) 



and taking the limit e — > 0, the time-derivative terms 



(100) 



and 
3 

W 



J dr^tr [ZU 



Or 



(101) 

result. These are identical to those derived microscopi- 
cally in Sec. HVBl 

We now discuss the phase diagram of the effective ac- 
tion. We shall find it convenient to work with discrete 
imaginary time (i.e. finite e). For simplicity, we shall 
set V c = t r c — t c and — t T s = t r s = t s , and we may set 
t s = 1 without further loss of generality, by rescaling the 
fermion fields. The phase diagram may now be drawn 
as a function of the dimensionless parameters t c and K\ 
this is shown schematically in Fig. [31 

We are most interested in understanding what happens 
when K is large. In this limit the gauge fluctuations are 
suppressed, and the spinons and rotor bosons become 
"good variables" with which to understand the physics. 
Furthermore, in this limit we will be able to access the 
SU(2) ASL, the semimetal, and a continuous quantum 
phase transition between them. 



K = < 



SU(2) ASL 


Mott transition 


Confined 
Insulator(s) 

(Neel, VBS, others...) 


Semimetal 



t = K=0 



t = oo 



FIG. 3: Schematic phase diagram of the effective low-energy 
theory S cS given in Eqs. j96j, (J9TJ) and (J98j) . The fermion 
hopping t s has been set to unity by rescaling the fermion 
fields. 



We first consider the limit K 3> 1 and ( c « 1. Here, 
the rotor bosons will be gapped, and we can safely in- 
tegrate them out to obtain an action involving only the 
fermions and the gauge field. At large K, the term Sg fa- 
vors gauge field configurations that are gauge-equivalent 
to U xx > = 1. If we set K to infinity we can drop the 
gauge field altogether in S c p , and the spinons behave as 
free Dirac fermions. For K large but finite, these Dirac 
fermions become coupled to a fluctuating SU(2) gauge 
field - we are therefore describing an SU(2) algebraic 
spin liquid. 

The simplest possibility is that the SU(2) ASL is a 
stable phase. While this assumption can be justified in 



the limit of a large number of fermion flavors;** we are 
not guaranteed that it holds for the Hubbard model. As 
we shall discuss in Sec. I VIII even if it is not stable, the 
SU(2) ASL may still have important implications for the 
metal-insulator transition. For the moment we simply 
assume that the SU(2) ASL is a stable phase. 

Next, we consider the limit K, t c ^> 1. The rotor 
bosons are condensed, resulting in a Higgs phase for 
the SU(2) gauge field. To understand the nature of this 
phase, we couple an external electromagnetic vector po- 
tential A xx i to the rotor bosons. Their part of the action 
then takes the form 



ooff _ 

— 



-tr 



(xx 1 ) 



tr e* 



Z\ UxT.' 'Z'X 



(102) 



is 



This form can be understood by noting that A 
minimally coupled to the r 3 -component of the pseu- 
dospin, which is simply the electromagnetic charge. Now 
we make the SU(2) gauge transformation (to "unitary 
gauge") 



U X x' 
Fr 



Z x U xx i Z^i 
F x Zl 



(103) 
(104) 
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The various terms in S now take the form 

S f = ~*c E to{e iA ™' r3 U xx ,) (105) 

{xx') 

Sf = J2 tofcU^F*,), (106) 

(xx') 

and the form of S c g s is unchanged. When t c is large, U xx > 
has only small fluctuations about e~ %A ^' T * , so we write 

U xx , = exp ( - U^r 3 ) exp ( ma: " 2 " r ) , (107) 

where d 1:E / may be treated as a small quantity. We can 
then expand Sf for small a xx i , and find (at leading non- 
vanishing order) a quadratic "mass" term for the fluctu- 
ations of the gauge field. Therefore the gauge field is 
massive in this phase, and can be integrated out. The 
fermion part of the action becomes 

Sf = ^2tr(F x e~ lA ^ T3 e lh ^'- T/2 F' x ) (108) 

(xx 1 ) 

« »(/•:,< ' " t3 K)- (io9) 

(xx') 

In the second line, the coupling to the massive gauge 
field has been dropped, as it will only generate irrele- 
vant four-fcrmion interactions. Therefore the low-energy 
theory consists simply of free Dirac fermions minimally 
coupled to the electromagnetic field. In fact, going to 
temporal gauge and taking the time-continuum limit we 
can rewrite this as 




which shows that these fermions are just electrons, and 
we are describing the semimetal. 

It should be noted that, even though we have con- 
densed the pseudospin-carrying field Z x , pseudospin is 
not broken in the Higgs phase, as is correct for the 
semimetal. To see this, we work for simplicity in the 
K = oo limit,— where U xx i is pure gauge, and we may 
choose U xx i = 1. Then, in the Higgs phase, 



( z xl z x2 \ 




Zl 


3 


V — Z x2 z xl J 









In this gauge, the "condensate" (Z x ) is independent of 
position x. By assumption, at least one of z~i,z 2 is 
nonzero, so Z x can be written as a real number times 
an SU(2) matrix. If instead we work in the gauge 
where U xx i = A X A^,, where are SU(2) matrices, then 
(Z x ) = A X Z X . Therefore, any two condensates can be 
rotated into one another by an appropriate SU(2) gauge 
transformation. 

Working in any fixed gauge, it might appear that the 
pseudospin symmetry has been broken, because Z x will 



transform nontrivially under pseudospin rotation. Con- 
sider a condensate specified by Z\ that can be trans- 
formed into Z x via a global pseudospin rotation. We can 
also transform Z\ into Z^. by making an SU(2) gauge 
transformation. Now, configurations that differ only by 
a gauge transformation should be thought of as different 
labels for the same physical stated therefore, Z x and 
Z x both describe the same state, and the apparent pseu- 
dospin breaking is a gauge artifact. 

Next, the Mott transition itself can be understood 
starting from the spin liquid phase. Barring a first-order 
transition, as t c is increased, the gap to rotor boson exci- 
tations will go to zero at a critical value of t c . The bosons 
condense on the other side of the critical point, and the 
semimetal is obtained. 

We also briefly discuss the other regions of the phase 
diagram, although these are not our main focus. For 
t c 3> 1 and arbitrary AT, we still expect to be in the 
semimetal phase - this is clearest if one works in unitary 
gauge as above. For AT, t c <C 1 the gauge field fluctuates 
strongly and the physics is less clear, especially given the 
coupling to the fermions. This region of the phase dia- 
gram is expected to be occupied primarily by insulating 
states where the gauge field is in a confining phase. The 
most likely such state, which is expected to occupy at 
least some part of the phase diagram, is the Neel antifer- 
romagnet. Valence-bond solids and other states are also 
possible. 

Finally, we are now in a position to comment on the 
analysis of Refs. I40ll4ll . where it was claimed that the 
SU(2) slave rotor formulation can lead to a description 
of a non-Fermi liquid metallic state. In terms appropri- 
ate to the present discussion and notation, that analysis 
identified the non-Fermi liquid as a state where, in some 
fixed gauge, (z\) ^ and (z 2 ) — 0. The Fermi liquid was 
identified as a state where both (z\) and (z 2 ) are nonzero. 
However, as per the discussion above, these two patterns 
of boson condensation can be rotated into one another by 
an SU(2) gauge transformation and thus do not describe 
distinct phases. This is true even away from half-filling 
where the pseudospin symmetr y is broken, which is the 
situation considered in Refs. \MM& In general, we feel it 
is probably not possible to describe a non-Fermi liquid 
conducting state in the SU(2) slave-rotor formulation. 



B. Field theory for the spin liquid 

We now construct a continuum quantum field theory 
for the SU(2) ASL; many results stated in this section are 
discussed in more detail in Appendix [Bl We begin with 
the continuum description of the spin liquid within mean- 
field theory. On the lattice the mean-field Hamiltonian 
is 

HMFT =t s J2 + H - C . ( 112 ) 

(rr f ) 
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where ip r is defined in Eq. (|40|) . We have included 
nearest-neighbor hopping only; as long as one allows ar- 
bitrary perturbations to the continuum effective action, 
this is sufficient to derive the low-energy theory. As dis- 
cussed in Sec. [I] and Fig- HI the SU(2) spinors tp r are also 
labeled as ijjRi , where R labels the unit cells of the hon- 
eycomb lattice and i = 1, 2 labels the two-site basis. The 
Fourier transform is then defined as 

where N c is the number of unit cells in the lattice. 

The Hamiltonian 7Ymft has gapless Dirac points at 
the Fermi level (fixed at zero energy by spinon particle- 
hole symmetry, which is part of the SU(2) gauge group). 
These nodal points lie at ±Q in the Brillouin zone, where 
Q = (iir/3)x. Working in momentum space, we define 
the continuum fields 

^ ~ ( ) < ii4 » 
M«)~(tz::)- (ii5 » 

where \q\ <C A, and A is a cutoff much smaller than the 
zone size. We can also Fourier transform to obtain real- 
space continuum fields 

Mr) = J J^S^M*), (H6) 

where the momentum cutoff A is implicit. 

The field tp a (where a = 1,2) is actually a four- 
component object, and it is convenient to define two 
kinds of Pauli matrices acting on it by defining the fol- 
lowing 4x4 matrices in block form: 

s (o °«) < 117 > 

T l = ( I 1 I 2 ) (118) 

V T 21 T 22 / 

V T 21P T 22P / 

We say that the fi l Pauli matrices act in the SU(2) gauge 
space, and the t 1 Pauli matrices act in the Dirac spinor 
or Lorentz space. 

The continuum mean-field Lagrangian density (in 
imaginary-time) is 

Co = <p a [ - ilftdf^tfa, (120) 

where (i = 0, 1, 2, we have defined 

£a=^T 3 , (121) 

and 

7/i = ( T 3, T 2,_ r i). (122) 



This is nothing but the Lagrangian for massless Dirac 
fcrmions in 2 + 1 dimensions. The coupling to the fluc- 
tuating SU(2) gauge field is introduced by replacing 
with the covariant derivative 

ia 1 ., u? 

D» = d lt + -f-. (123) 

We also define the field-strength tensor for the gauge field 
as 

fU = dA - d v ai + e«*a£a*. (124) 
The effective Lagrangian for the SU(2) ASL is then 

a ff SL = M- hm^K + ^fUU- ( 125 ) 

The gauge coupling constant g 2 has dimensions of mass. 
Other terms consistent with the symmetries of the under- 
lying model are added as perturbations, and local stabil- 
ity against such perturbations can be considered. 

As it stands this field theory is strongly coupled and 
cannot be analyzed reliably - a straightforward per- 
turbation theory leads to severe infrared divergences. 
However a controlled analysis is possible in a large-A/ 
expansion ) 35 i 38 where instead of a — 1,2, one allows 
a = 1 , . . . , Nf . (In our notation, the physical case corre- 
sponds to Nf = 2.) The perturbation theory is reorga- 
nized into an expansion in powers of 1/Nf, whereupon 
the infrared divergences are cured and systematic cal- 
culation becomes possible. As in the case of the U(l) 
ASL ; 26 ' 27 ' 34 one arrives at the conclusion that the SU(2) 
ASL is an interacting critical state, with nontrivial criti- 
cal exponents that vary as a function of 1/Nf. Also as in 
the U(l) case, 2 - this state is a stable phase in the large- 
Nf limit; the only potentially dangerous perturbations 
are fermion mass terms or "velocity anisotropy" terms 
(i.e. fermion bilincars with one derivative), and on the 
honeycomb lattice it is easy to show that all such terms 
are forbidden by symmetry. 

Very recently, in the context of an SU(2) ASL in the 
square lattice Heisenberg model (the 7r-flux stat o 24 ' 25 ), 
Ran and Wen have shown that the SU(2) ASL has an 
emergent Sp(2Ay) symmetry? 3 ^ This symmetry is present 
exactly in Eq. (|125[) . but it is hidden and some work is 
needed to expose it (see Ref. HH and Appendix [B] for de- 
tails). The Sp(2iV/) plays a role analogous to the SU (Nf) 
emergent symmetry recently discussed in the context of 
U(l) algebraic spin liquids-^ in particular, it also leads 
to a unification of seemingly unrelated "competing or- 
ders." These are a set of observable quantities that have 
slowly decaying power-law correlations, all with the same 
critical exponent. Ref. [H found that a set of five fermion 
bilinears are unified together into an SO (5) vector repre- 
sentation of the Sp(4) symmetry (recall that Sp(4) and 
SO(5) have the same Lie algebra), and that coupling to 
the SU(2) gauge field results in slowly decaying correla- 
tions for this SO (5) vector. For the square lattice Heisen- 
berg model, they found that three components of the 
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FIG. 4: Depiction of a valence-bond solid state on the hon- 
eycomb lattice. The shaded regions represent spin-singlet 
valence bonds. This VBS pattern corresponds to one of 
the fluctuating competing orders present in the SU(2) ASL; 
the other such competing order is simply Neel antiferromag- 
netism. Both these observables have slowly decaying power- 
law correlations characterized by the same critical exponent; 
this phenomenon should be a useful test for the SU(2) ASL 
in future numerical studies. 



the critical exponents for Neel and VBS correlations are 
the same, the simultaneous presence of slowly decaying 
correlations in both observables would already be a sign 
of the SU(2) ASL. 

It may be useful to mention that VBS correlations 
can be measured in two different ways. The first is, of 
course, to look at a dimer-dimer correlation function as 
in Eq. (|128p . However, because there is no difference at 
the level of symmetry between the two quantities, one 
can just as well look as correlations of the bond kinetic 
energy density. This would be simpler in determinantal 
quantum Monte Carlo simulations, because the relevant 
correlation function involves only four electron operators, 
as opposed to eight in the dimer-dimer correlator. The 
kinetic energy on the bond (r, r') is 



K r 



- -(A 



(<Z a c r , a + H.c), 



(129) 



and an example of its long-distance correlation function 
in the SU(2) ASL is 



(K(m,R2)K( 01t 02)) - ((-^(01,02))) 



2 cos(4tt.R x /3) 



\R 



2A 



with the same critical exponent A as above. 



(130) 



vector correspond to the Neel order parameter, and the 
other two to the order parameter for columnar and box 
valence-bond solid (VBS) states^ 

The result for the honeycomb lattice is very similar; 
the SO (5) vector again contains the Neel order parame- 
ter n and the two-component order parameter for a VBS 
state, v — (v\ : V2)- In Appendix IBl these fields are de- 
fined as gauge-invariant bilinears of the ip a fermion fields. 
The VBS order parameter corresponds to the pattern of 
dimerization shown in Fig. 0] Using the result of Ref. [H, 
to leading order in the 1/Nf expansion the scaling di- 
mension of these operators is 



A = dimjn] = dim^] = 2 



32 
n 2 N f ' 



(126) 



This means that Neel correlations, and the VBS correla- 
tions pictured in Fig. [?J fall off with the same power-law 
decay at long distances. For instance, 



>R1 



Soi 



1 



\R 



2A ' 



(127) 



and 



2 cos(47ri? x /3) 



l-Ri 



2A 



({Sri ■ Sr2){Sqi ■ S02)) — {{S01 ■ S02)) 

(128) 

It should be noted that the magnitude of the 1 /Nf cor- 
rection to A is rather large, so the above correlators may 
decay quite slowly, making them easy to observe in nu- 
merical simulations. Furthermore, even if the system 
sizes available are not sufficient to determine whether 



C. Field theory for the Mott transition 

To construct the continuum field theory for the Mott 
transition, we start with that discussed above for the 
spin liquid state, and add in the rotor bosons. To under- 
stand the form of the bosonic sector, we again start at 
the mean-field level. On the lattice, the bosons obey the 
mean-field action 



J MFT 



5^(tr(4 +ez ^ 



{ tl i Z lr,r) Z (r',T) 

1 +' ■"^wz,), 



r.r 



(131) 



where for simplicity the coefficients of the first two terms 
have been set to unity, and 



Zxl 



Zx2 
z x2 z xl 



(132) 



where z x i, z X 2 are arbitrary complex numbers. The con- 
stant r c o is chosen so that the lowest energy mode of the 
boson field just becomes gapless as r is decreased to zero. 
It is a straightforward task to diagonalize the quadratic 
form of Eq. p31[) . and one finds that for r > there is 
a single low-energy mode at k = in the Brillouin zone. 
The low-energy fluctuations can be described by a con- 
tinuum 2x2 matrix field O, which can be related to the 
lattice fields by 



6(ir) ~ Z[r,(H,l)] + Z [t,(R,2)]- 



(133) 
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Here, (R, 1) and (R, 2) label the two honeycomb sites in 
the unit cell at R. We write 



e 



where 9i, Qi are arbitrary complex numbers. 
The mean-field Lagrangian density for O is 

£? = ^tr (e+a> M e) + ^tr(ete), 



(134) 



(135) 



where, in the first term, the first partial derivative acts 
to the left. We have set the velocity of the O boson to 
unity - we come back to this point shortly. It should be 
noted that no linear time derivative term is present. Such 
a term is forbidden by particle-hole symmetry (together 
with invariance under global SU(2) gauge transforma- 
tions and U(l) charge rotations). This has the important 
consequence that the bosonic sector obeys relativistic dy- 
namics at low energy. 

The coupling to the gauge field is included by intro- 
ducing covariant derivatives, and the kinetic energy term 
for 9 becomes 



^ 1 / x ,<— id 1 ,, a 1 . , 



ia l .,a l , 



e . (136) 



The full effective action describing the Mott transition is 
then 



£° M ff o tt = AIl + £L + \ tr (ete) + \ ( tr (©to) 

(137) 

where we have also included a quartic term for the 
boson field. This can be written in a more conventional 
form if we define the two-component field 

= ( °L V (138) 



The bosonic part of the Lagrangian then becomes 

£g n + ^tr(ete) + ^(tr(0te); 2 



(9, 



)tf + n? f tf + A(tf f t?) 2 (139) 



So the full field theory is simply Nf — 2 Dirac fermions 
and Nb = 1 relativistic bosons, both coupled to an SU(2) 
gauge field. All the matter fields transform in the funda- 
mental representation of the gauge group. 

The Mott transition from the spin liquid will occur 
when r is tuned through a critical value r c and the mass 
for vanishes. On the conducting side of the transition, 
condenses and gives the gauge boson a mass via the 
Higgs mechanism. As discussed above in Sec. IVI A[ this 
results in a description of the semimetal. The critical 
properties of the Mott transition can be analyzed in at 
least two controlled limits: one can make the number 
of both fermion (Nf) and boson (Nb) fields large and 



perform a simultaneous expansion in 1/Nf and 1/Nb, or 
one can carry out an expansion in 4 - e space-time di- 
mensions. As our primary focus is on signatures of the 
spin liquid phase itself, we have left such calculations for 
future work. If evidence for spin liquid physics of the 
type we propose here is indeed found in the honeycomb 
Hubbard model, it would be important to work out the 
properties of the critical point in more detail. 

A word about Lorentz invariance is now in order. 
There are three bare velocities in the effective Lagrangian 

- that of the fermions, the boson and the gauge boson 

- and in order to have Lorentz invariance they must all be 
equal. The bare gauge boson velocity is actually not im- 
portant, because the long-distance, low-energy dynamics 
for the gauge bosons are generated entirely by the matter 
fields. One of the remaining velocities can be set to unity 
by an appropriate choice of units, but the ratio of the 
fermion and boson velocities is a dimensionless parame- 
ter that cannot be removed by scaling, and in principle 
may be important at the Mott critical point. Supposing 
the difference in velocities is small, it may be treated as 
a perturbation starting from the Lorentz-invariant criti- 
cal point, which could be irrelevant, marginal or relevant 
in the RG sense. While this can be resolved by RG cal- 
culations of the type described above, the answer is not 
immediately obvious, because, in both the Nf,Nb — > oo 
and e — > limits, the difference in velocities is marginal, 
and its fate will be determined by perturbative correc- 
tions in the appropriate small parameters. The most 
likely outcome, which occurs in other caseS) 36 i 45 is that 
the difference in velocities will be irrelevant and the sys- 
tem will flow to the Lorentz-invariant critical point. This 
is an important issue; if the difference in velocities is rele- 
vant, then even in the large- N or small-e limit the critical 
point will have two relevant directions and will be unsta- 
ble. While we feel this is a very unlikely scenario, this 
issue will need to be addressed in any future study of the 
critical properties. 



D. Pseudospin breaking and U(l) ASL 

The SU(2) slave-rotor formulation has allowed us to 
understand the Mott transition out of the semimetal into 
the SU(2) ASL, so far assuming the presence of pseu- 
dospin symmetry On the other hand, the U(l) slave- 
rotor formulation finds instead a Mott transition into a 
U(l) ASLJ£ While the U(l) formulation does not respect 
the pseudospin symmetry, we might expect to recover its 
results in the presence of perturbations that break pseu- 
dospin down to the 0(2) subgroup generated by U(l) 
charge rotations and particle-hole symmetry. It is inter- 
esting to explore this possibility because the U(l) ASL is 
likely to be more stable than its SU(2) counterpart, and 
also because it possesses a different pattern of compet- 
ing orders with slowly decaying power-law correlations, 30 
and can thus be easily distinguished from the SU(2) ASL. 
Furthermore, this situation should be accessible to quan- 
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turn Monte Carlo simulation without a sign problem, for 
example by adding the nearest-neighbor interaction "Hy 
ofEq. 

Here, we consider the effect of adding a pseudospin- 
breaking term to the low-energy theory of Sec. IVI Al 
The term we choose to add is not simply Hy written 
in slave-rotor variables; however, its symmetry transfor- 
mation properties are identical to those of Hy, and it will 
be generated upon integrating out some high energy de- 
grees of freedom. Other terms will of course also be gen- 
erated, so we cannot say for sure whether Hy will have 
the effects discussed here. Starting from the SU(2) ASL 
phase, we find that when the pseudospin breaking is large 
enough a quantum phase transition to the U(l) ASL re- 
sults. This phase transition is interesting in its own right 
as an example of a critical point with no order parame- 
ter on either side of the transition,— and, while we have 
not pursued this here, it may be fruitful to investigate the 
critical properties in more detail. The U(l) ASL can itself 
undergo an insulator-metal transition (to the semimetal 
state) upon condensation of a single bosonic field carry- 
ing unit charge under the U(l) gauge field. In the sector 
with strong pseudospin breaking, we therefore recover the 
effective theory obtained in the U(l) formulation^ 

It will now be convenient to introduce some additional 
notation. We define the composite field 

[Z(x)).. = \tT{a i Zt^Z x ). (140) 

It can be shown that Z(x) is an SO(3) matrix, and that, 
furthermore, any SO (3) matrix can be expressed in this 
form. This mapping is two-to-one from Z to Z (because 
Z and —Z map to the same image); in fact, Eq. (|140p 
is nothing but the well-known "double cover" map from 
SU(2) to SO (3). Under an SU(2) gauge transformation 
sending Z x — > e~ ta,c '' T ' 2 Z X) we have 

Z(x) -Z(as) [R(a x )] T , (141) 

where R(a x ) is the SO (3) matrix describing a rotation 
about the a^/la^l axis by an angle \a x \. Under a pseu- 
dospin rotation Z x — > Z x e la rT / 2 , Z undergoes a left- 
SO(3) rotation: 

Z(x) -> R{a)Z{x). (142) 

Next, we make use of the double-cover again to define 
an SO (3) gauge field: 

[<3 X % = \*{VL><r i U**>o*). (143) 
Under the SU(2) gauge transformation U xx > — > 

e -io B -«r/2 E ; a;a! , e »a a ,/- -/2 ) we haye 

G xx ' -> R(a x )Q xx ' [R(a x ,)] T . (144) 

The desired pseudospin-breaking term is a gauge- 
invariant kinetic energy for the composite Z boson: 

Sf = -td tr (R{2A xx ,z)M Z(x)Q xx ' [Z{x')] T ). 

{xx') 

(145) 



Here, M is a symmetric 3x3 matrix, and we have in- 
cluded the coupling to the electromagnetic vector poten- 
tial A xx i via the rotation matrix R(2A xx iz). If we choose 
M. = 1, then this kinetic term is pseudospin invariant. 
Instead, we choose M. = diag(0, 0,1), which breaks the 
pseudospin down to 0(2) as discussed above. For this 
choice of M, the dependence on A xx i drops out, and we 
are free to set A xx > = in Eq. (|145[) . So the only depen- 
dence on A xx , is in Sf [see Eq. (fT02|) ]. 

Now we consider the effect of making t c % ^> 1, making 
no assumptions for the moment about the other param- 
eters t c and K. We may proceed by expanding in small 
fluctuations about the minimum of S% ■ . We make the 
following transformation: 

U m , -> Z X U XX ,Z\, (146) 
F x - F x Zt (147) 

which also implies 

g xx ' -> [z( x )] T g xx 'z(x'). (148) 

We then have 

Sf = -t c2 tr (M Q xx ') ■ (149) 

{xx') 

Choosing Q xx (and hence U xx i) to minimize Sf > we find 
that 

U xx > = exp (ia xx >cr 3 ). (150) 

We have thus obtained a Higgs phase in which the 
SU(2) gauge field has been broken down to U(l). This 
corresponds to a condensation of certain components of 
the composite object Z(x). To better understand the 
nature of this phase, we work with the remaining terms 
in the action. We have 

Sf = -2te C0S ( A **' + a **') ( 151 ) 
{xx') 

Sf = ^ ir (/;,."•• ' 7 V:.) (152) 

{xx') 

Sf = -4FT^ cos ([V x a] p ). (153) 
p 

Here, [V x a] p is the lattice curl of a xx i, taken around 
the perimeter of the plaquette p. Next, we make the 
transformation a xx i — * a xx r + 6 X — x i and also F x ~> 
F x e~ l0 * a ; this introduces the bosonic matter field 9 X 
and restores a U(l) gauge invariance. The form of Sf 
and Sf is unchanged, and we have 

Sf = -2t c £ cos (6 X - 9 X , + A xx , + a xx ,). (154) 

{xx') 

We have thus arrived at the effective action obtained 
in the U(l) slave-rotor formulation. The compact 
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FIG. 5: Schematic phase diagrams of the effective action in 
the t c -t C 2 plane; K is taken to be large. As discussed in 
the text, two scenarios are shown, with differing assumptions 
about the critical behavior at the SU(2) ASL to semimetal 
transition, and hence different topologies of the phase dia- 
gram. 



U(l) gauge field a xx > is coupled to the charge- neutral, 
fermionic spinons and also the rotor boson field with 
phase variable 9 X . The 0-boson carries unit electromag- 
netic charge, as is apparent from the presence of the ex- 
ternal electromagnetic field in Sg. 

We now consider the consequences of the above result 
for the phase diagram, focusing on the regime K 1. 
It is useful to imagine starting in the SU(2) ASL phase, 
where t Cl t C 2 <C 1. As t C 2 is increased, eventually Z will 
condense and the U(l) gauge theory action given above 
is appropriate. Because t c <C 1, the #-boson is gapped, 
so the resulting phase is the U(l) ASL. The quantum 
critical point where Z condenses is thus a phase transi- 
tion between the SU(2) and U(l) algebraic spin liquids. 
Next, we keep t C 2 large and increase t c . Eventually the 
#-boson will condense, leading to the semimetal^ So, in 
the regime t C 2 ^ 1, we recover the scenario of Ref. Ha- 
Also, it should be noted that the semimetal phase per- 
sists into the regime t c ^> 1 and ( C 2 « 1. 

Schematic phase diagrams reflecting these considera- 
tions are shown in Fig. [5j Depending on the (unknown) 
critical properties of the SU(2) ASL-scmimctal transi- 
tion, different topologies of the phase diagram are pos- 
sible; we discuss the two simplest scenarios. In Fig. [5jV, 
pseudospin breaking is assumed to be relevant; in this 
scenario, the Mott transition for t C 2 ^ is always the 
U(l) ASL-semimetal transition. On the other hand, in 
Fig. [SB, pseudospin breaking is taken to be irrelevant, 
and t C 2 must be increased above a critical value to change 
the nature of the Mott transition. Other, more compli- 
cated scenarios are of course also possible. It should be 
noted that, starting from within the SU(2) ASL phase, 



t C 2 must always be increased above a critical value to 
reach the phase transition to the U(l) ASL; this follows 
from the Z-boson gap in the SU(2) ASL. 



VII. DISCUSSION 

The SU(2) slave-rotor formulation developed in this 
paper is a general formalism to describe spin liquids and 
related physics in Hubbard models; it is especially useful 
at half-filling to describe spin liquid ground states near 
the Mott transition. We proposed the honeycomb Hub- 
bard model as a candidate for this physics, and showed 
that the simplest possibility is a direct, continuous quan- 
tum phase transition between the conducting semimetal 
and an SU(2) algebraic spin liquid. 

We hope that this physics will be looked for in fu- 
ture quantum Monte Carlo simulations of the honeycomb 
model; to that end we summarize our main predictions 
and suggestions for such work. First, and most simply, 
the mean-field analysis suggests that the spin liquid phase 
is more likely to be found if a third-neighbor electron 
hopping t" is present, when t"/t < 0. Second, a robust 
signature of the SU(2) ASL is crucial for its detection 
in numerical simulations, especially as one may be lim- 
ited to rather small systems, and only a small region 
of spin liquid phase may be present. The simultaneous 
presence of slowly decaying power-law Neel and valence- 
bond solid correlations, with the same critical exponent, 
should provide such a signature. The particular pattern 
of VBS order with slowly decaying correlations is shown 
in Fig. [31 and some specific examples of correlation func- 
tions exhibiting the slow power-law decay are given in 
Sec. IVIB1 The VBS correlations can be measured ei- 
ther using a dimer-dimer correlator, or by the correlation 
function of the bond kinetic energy. While limitation to 
small system sizes may prevent a detailed analysis of the 
critical behavior, it should be possible to look for the 
simultaneous presence of slowly-decaying Neel and VBS 
correlations, which would already be some evidence for 
the SU(2) ASL. 

Finally, it may be useful to consider an extended Hub- 
bard model with the nearest-neighbor density interaction 
Tiv of Eq. ((2]). This term breaks the pseudospin symme- 
try down to 0(2), and, as discussed in Sec. IVIDl it may 
favor a U(l) ASL over the SU(2) ASL. This is poten- 
tially advantageous because the U(l) ASL is likely more 
stable against small perturbations than its SU(2) coun- 
terpart, because gauge theories with more gauge bosons 
(i.e. more colors) have a stronger tendency toward con- 
finement. 

Throughout this paper we have been assuming that the 
SU(2) ASL is a stable phase; while this does hold in the 
large- Nf limit, it may not be true in the physical case 
of the Hubbard model. One possibility is that the SU(2) 
ASL is not a stable phase, but its U(l) counterpart is, and 
in that case a pseudospin breaking term such as TLv may 
be required to obtain a stable spin liquid adjacent to the 
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Mott transition in the honeycomb Hubbard model. (We 
caution that the U(l) ASL itself is also not guaranteed 
to be a stable phase.) 

As in the SU(2) case, the U(l) ASL can be tested 
for by looking for competing orders with slowly decay- 
ing power-law correlations. These competing orders are 
unified together by an SU(4) symmetry, and their cor- 
relations fall off with the same critical exponent - the 
value of this exponent will be different from its counter- 
part in the SU(2) ASL. The pattern of competing orders 
could be determined in a straightforward fashion follow- 
ing Ref.Hil and is a superset of those found in the SU(2) 
ASL, including the Neel and VBS correlations discussed 
here, as well as correlations in other observables. 

For the case of nearest-neighbor hopping only, the 
honeycomb Hubbard model has actually already been 
studied via quantum Monte Carlo and series expansion 
techniques. 31 Those results suggest, but do not conclu- 
sively establish, a direct, continuous quantum phase tran- 
sition between the SM and AF states for (U/t) c ~ 4 — 5, 
which could be the conventional SM-AF transition 32 
shown in Fig. [2^. Of course, a small region of spin liquid 
intervening between SM and AF is not excluded. An- 
other possibility, if the SU(2) ASL is not a stable phase, 
is a deconfined critical poin t 46 i 47 from SM to AF. Briefly, 
this can be understood as a transition from SM to the 
SU(2) ASL, which is in turn unstable (via a dangerously 
irrelevant perturbation) to Neel order. In order to obtain 
some information about which scenario is in fact occur- 
ring, the dimer-dimer or bond kinetic energy correlation 
functions could be measured to look for an enhanced ten- 
dency toward VBS ordering near the transition. Also, 
as discussed above, adding a t" hopping or the nearest- 
neighbor repulsion Tly may help to stabilize a spin liquid 
phase or deconfined critical point. 

While our main focus in this paper has been the hon- 
eycomb lattice, the SU(2) slave-rotor formulation should 
be useful much more generally, as an extension of the 
SU(2) gauge theory of the Heisenberg model to the case 
of the Hubbard model. It would clearly be interesting to 
investigate the triangular lattice using this approach, mo- 
tivated by experiments on k-(BEDT-TTF) 2 Cu 2 (CN) 3 . 
This would allow the competition among a wider range 
of spin liquids than that considered in Ref. Ojl to be ad- 
dressed. The simplest way to carry out such an analysis 
is via mean-field theory, but a study of variational wave- 
functions is probably more reliable. It may be possible 
to use the SU(2) slave-rotor formulation to develop new 
variational wavefunctions for Hubbard models, to study 
physics near the Mott transition. 

On an even more general note, the SU(2) slave-rotor 
formulation may help to shed more light on the physical 
nature of fermionic spinons. For most spin liquids with 
fermionic spinons, we lack a simple physical picture of 
the spinons in terms of microscopic degrees of freedom. 
One can draw the usual cartoon pictures of resonating 
valence-bond states, and show that breaking a valence 
bond corresponds to creating a pair of spinons; however, 



such pictures do nothing to distinguish among the great 
variety of different spin liquid states, and do not even tell 
us whether the spinons should be thought of as fermions 
or bosons. 

It turns out that more progress has been made in com- 
ing to a physical understanding of bosonic spinons. The 
reason for this is that spin liquids with bosonic spinons 
can be thought of as adjacent to a magnetically ordered 
state, obtained upon condensation of the spinons. In 
some cases, the spin liquid phase, and hence the spinons, 
can then be understood in terms of the degrees of freedom 
of the ordered state. For example, in two-dimensional 
5 = 1/2 systems with XY symmetry, a Z 2 spin liquid 
be obtained upon condensation of pairs of vortices in 
the XY-ordered stated The deconfined critical point be- 
tween the VBS and the XY-ordered state, and its bosonic 
spinons, can also be understood in terms of the vortices 
of the XY-ordered state^&il 

The SU(2) slave-rotor formulation suggests a similar 
viewpoint for spin liquids with fermionic spinons, as it 
should allow for a description of the Mott transition out 
of any such spin liquid, into a conventional conducting 
state (either a metal or superconductor). At the tran- 
sition, the fermionic spinons transform into electron-like 
quasiparticles. This suggests that a physical picture of 
the spinons might be obtained by understanding the spin 
liquid in terms of the degrees of freedom of the nearby 
conducting state. This may be a challenging task, but it 
has been carried out in one case: starting from a d-wave 
superconductor with gapless nodal quasiparticles, con- 
densation of pairs of superconducting vortices leads to a 
Zi spin liquid with nodal fermionic spinons i 48 ' 49 i 50 ' 51 We 
are optimistic that further progress along similar lines is 
possible in other cases. 
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APPENDIX A: QUANTUM MECHANICS OF 
THE SU(2) MATRIX ROTOR 

In this appendix we shall derive some details of the 
quantum mechanics of the SU(2) matrix rotor. In par- 
ticular, we develop its functional integral representation, 
both with and without an SU(2) gauge constraint. As a 
warmup we first review the 0(3) quantum rotor, then, 
in the following section we return to the SU(2) matrix 
rotor. 
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1. 0(3) quantum rotor 

Consider a particle constrained to move on the surface 
of the 2-sphere S 2 . Its position is given by the unit- 
length 3-vector x (i.e. x 2 = 1). The quantum Hilbcrt 
space has a basis of position eigenstates {|a;}}, with the 
usual inner product (x\x') = 5(x — x'). Note that this is 
a "two-dimensional" delta function, as appropriate for a 
variable constrained to a 2-dimensional manifold, so, for 
example, j g2 dx 5(x) = 1. 

Angular momentum operators L are defined as the 
generators of infinitesimal rotations. That is, for a finite 
rotation, we have 



x) = \R(a)x), 



(Al) 



where R(a) denotes a rotation by |a| about the a/|a| 
axis. This implies the familiar commutation relations 
[L\U\ = ie^ k L k and [L\x j ] = ie ijk x k . 

Next, we construct a basis of angular momentum eigen- 
states. First of all, there is a unique singlet state given 
by 



|0> = -J= / dx\x). 
V47T Js 2 



(A2) 



The higher angular momentum eigenstates are spherical 
harmonics, which can be obtained by acting on |0) with 
irreducible polynomials of xf, that is, polynomials that 
are symmetric in 0(3) vector indices and traceless in any 
pair of two such indices. Such polynomials transform 
under irreducible representations of the rotation group, 
hence the name. For example, we have, not worrying 
about normalization, 



\i) = Xi \0) 

\ij) = {xiXj - -%)|0) 



(A3) 
(A4) 



\ijk) = [xiXjX k - -(SijXk + 8i k Xj + SjkXi)] |0)(A5) 

These states transform with total angular momentum I = 
1,2 and 3, respectively. Including arbitrary high values 
of £, the angular momentum eigenstates form a basis for 
the Hilbcrt space. 

Consider the Hamiltonian 



H=-L 2 . 

2 



(A6) 



We would like to go to a functional integral description of 
the partition function. In order to do this, we will need 
to evaluate objects like 



(x(r + e)\e-^\x(r)) 



(A7) 



As usual we assume only smooth configurations of x(t) 
contribute to the path integral, so that x(r + e) = x(t) + 
0(e). It should be enough to consider only the first- 
order term in e, which will dominate the 0{e 2 ) and higher 



terms. Then, to first order in e, we note that x(t + e) 
lies in the plane tangent to the sphere at x(t). So, for 
the purposes of evaluating this matrix element, we can 
change the Hilbert space to that of a particle moving in 
this tangent plane. So we replace x — > xt, where Xt lies 
in the plane tangent to the sphere at x(t), and Xt{t) = 
x(t). We then have 

(x(r + e)\e- eH \x(r)) - (x t (t + e)\e^ HT \x t {t)). (A8) 

In order to complete this we have to understand how to 
define Ht, the Hamiltonian in the tangent-plane Hilbcrt 
space. This is seen most easily by considering a concrete 
example, where x(t) = z. Then we see that 



JBL* 



e ieLX \z) 



i6L y 



z) = \z) (A9) 
\z-dy + 0(d 2 )) (A10) 
\z) = \z + 6x + 0(6 2 )). (All) 



For infinitesimal rotations, the angular momentum oper- 
ators generate translations in the tangent plane, and so 
should be replaced with the linear momentum operators 
in the tangent plane Hilbcrt space. Specifically we have 



L y 



-Py 



Px 

0. 



(A12) 
(A13) 
(A14) 



More generally, if we wish to consider the tangent plane 
at an arbitrary point x , we have 



Xq X p_L, 



(A15) 



where p± is the linear momentum in the tangent plane, 
and p± ■ Xo = 0. Under this replacement, note that we 
have simply L 2 — > p 2 ± , so we have H T = \p\ ■ 



Thus we can write 



(x(t + e)\e- eH \x(r)) « (x T (r + e)\e'^\x T ( T )) 



1 



d 2 p ± (x T (T + e)\ Pl _}e-^ p ^(p^\x T (T)). (A16) 



It is now clear that we will arrive at the functional inte- 
gral: 



J Vx{t) j V P± (t 7 x{t)) ( 



-S[x,p ± ] 



(A17) 

where x is constrained to lie on the unit sphere, it should 
be noted that the range of integration of p± depends on 
x(t), and the action has the form 

S[x,p±]= J dr(- ip ± ■ x + ^pI). (A18) 

Upon integrating out the momentum variable we obtain 
the familiar form 

S[x] = lj J o dr(d T x) 2 . (A19) 
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2. SU(2) matrix rotor 

Now that we have dispensed with our warmup, let us 
consider the situation we need for the analysis in the pa- 
per, the quantum mechanics of the SU(2) matrix rotor. 
As mentioned in Sec. IIVB1 this is identical to a parti- 
cle constrained to the 3-sphere with 0(4) symmetry, and 
here we shall find it convenient to work with the 0(4) 
rotor description, especially to derive the functional in- 
tegral. The position of the particle is specified by the 
SU(2) matrix Z, or equivalcntly by the unit-length 4- 
vector X/j,, with /i = 0, . . . , 3. The equivalence of these 
two representations if manifest upon writing 

Z = x - icr ■ x, (A20) 

where x = (x\, X2,x^). The Hilbert space is specified by 
the basis of position eigenstates or, equivalcntly, 

{!*>}• 

SO (4) rotations of the position vector are gener- 
ated by the operators = —L vti , which satisfy the 
commutation relations 

[L^vi L\ a ] = i(5fj,\L ucr + 5 va L^\ — S^ a L u \ — S u \L^ a ), 

(A21) 

and 

[Ln U ,x\] = iiSpxXu - 8 v xx^). (A22) 

Now, SO(4) = SU(2) x SU(2)/Z 2 , and we can express the 
generators of right and left SU(2) rotations of Z, Jr and 
Jl, in terms of the L^. To see this, we follow the usual 
procedure and define 

J* = 2 e ijkLjk (A23) 
Ki = L 0l , (A24) 

where i = 1, . . . , 3. In terms of these operators 

Jr = \( J - K ) ( A25 ) 
Jl = \{J + K), (A26) 

which satisfy the commutation relations and other iden- 
tities given in Sec. IIV Bl 

As for the 0(3) rotor, we need to define a basis of 
angular momentum eigenstates. There is again a unique 
rotationally-invariant state 

|0) rot = -2= / d 3 x\x). (A27) 
V^T Js a 

By analogy with the 0(3) rotor, we should be able to 
access all higher angular momentum states by acting on 
|0) with irreducible polynomials in x^. Again, such poly- 
nomials are completely symmetric in 0(4) vector indices, 
traceless over any pair of indices, and therefore transform 
as irreducible representations of 0(4). For example, x^ 



and [Xfj,x u — (1/4)5^] satisfy these conditions. For the 
irreducible polynomial of degree k, there are dim^ inde- 
pendent components. A simple counting argument shows 

dim fe = (fc + l) 2 . (A28) 

Because they transform as irreducible representations 
of the S0(4) and hence SU(2) x SU(2) Lie algebra, the 
states just constructed can be labeled by the quantum 
numbers (£r,£l), the total angular momentum quan- 
tum numbers under right and left rotations, respectively. 
Only some pairs (£r,£l) actually appear in the Hilbert 
space - we will show that the multiplets with £r = £l 
form a complete basis. Now, we observe that we can 
construct the representation (£, £) in our Hilbert space, 
for £ — 0,1/2,1,3/2,.... This is done most easily in 
terms of which we can write as a two-index tensor 
Z aa , where I use roman letters for the first index and 
Greek letters for the second one because they transform 
under independent SU(2)'s. The state corresponding to 
(0, 0) is just |0) ro t, and clearly the (1/2, 1/2) multiplet is 
made up of the states Z aa \0} rot . In general, forgetting 
about normalization, we have 

(n/2, n/2) 4-> ]T Z„ p(ljai • • • Z ap(n)0in |0) rot , (A29) 
p 

where the sum is over all permutations of {l,...,n}. 
These states are completely symmetric in both the ai 
and Ui indices, which implies that they transform as a 
(n/2, n/2) multiplet. Now, the (n/2, n/2) multiplet con- 
tains (n + l) 2 states, so we are forced to conclude that 
the states of Eq. (|A29|) correspond to those obtained from 
the irreducible polynomial of degree n. This exhausts all 
possible states, so the (£, £) multiplets form a basis. 

Next we will construct a functional integral, first for 
the Hamiltonian 

H= I -^2 L % = I (J 2 R + Jl)- (A30) 

Again, to evaluate each matrix element in the Trotter 
expansion, we need to replace the Hilbert space of the 
3-sphere by the Hilbert space of the tangent hyperspace. 
So, for a particle with position x, we make the replace- 
ment L^v — > x^p v — x v p^, where p^x^ = 0. Under this 
replacement we have Y) <v L 2 „ — > p^p^. It is then clear 
that the functional integral takes the same form as be- 
fore, but in one higher dimension: 

Z = J Vx{t) J Vp(T,x(T))e- s[x ' p] (A31) 

S[x,p] = J dr^-ip^x^+^p^p^. (A32) 

Again, upon integrating out the momentum variable, we 
have 

1 r p 
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Let us rewrite this in SU(2) matrix notation. We have 
tx{d T Z^d T Z) 



tr 



(d T Xo + icr ■ d T x)(d T xo — icr ■ d T x) 



= 2(d T x^) 2 , 
and therefore 



VZ(r)e 



-S[Z] 



i r 

S[Z] = — / dr tr(d T Z^d r Z). 



(A34) 



(A35) 



(A36) 



Here, Z(t) is constrained to be an SU(2) matrix. 

Now let us consider a different Hamiltonian problem, 
which is exactly what we will need to write down the 
functional integral form of the Hubbard model in the 
SU(2) slave rotor formulation. Suppose that H = IJ%, 
and there is a constraint involving J^: 



Jl + v = 0. 



(A37) 



Here, v is a vector that may depend on other dynamical 
fields (e.g. the fermionic spinons) but commutes with all 
operators acting in the rotor Hilbert space. So for the 
present considerations we can take it to be a c-number. 
We will define a functional integral by tracing over the 
unconstrained Hilbert space, but we define a projection 
operator V that projects onto the constrained Hilbert 
space. Then we have 



Z = tr{Ve- pH ) = tT(Ve~ e "Ve 



eH-n-eH 



Ve- eH ), (A38) 



which holds because [P,H] — 0. Now we have 
V cx / d\e~ leX<jL+v \ 



(A39) 



where the unimportant constant of proportionality will 
be absorbed into the functional integration measure. We 
now see that we need to evaluate objects like 



(x(t + e)|e 



e\-(J L +v) -eH 



\x(t)), (A40) 



which can be evaluated precisely as above by going to 
the tangent hyperplane Hilbert space. It can be shown 
that the correct replacement for J\ is 

1 



J R -> iPuPl*- 



(A41) 
We then obtain 

Z = J Vx(t)VX(t) J Vp(T,x{T))e- s ^ x ' p > x \ (A42) 
and 



S[x,p,X] 



dr(- ip^Xf, + jPfiPp. 



(A43) 



-A • (x x p + x p — xp ) + iX ■ v 



Upon integrating out the momentum, we obtain 



S[x,X] 



13 



dr 



j{d T x + ■ x) 



1 



1 



— (d T x — -A xx— —xqX) + iX ■ v 



21 



dri 



{tr \z^(d T 



iX ■ cr 



){dr 



(A44) 
^)Z] 



+ iX ■ v\, 



(A45) 



where the second equality can be verified by simple alge- 
bra after expressing Z = xo — icr ■ x. 



APPENDIX B: SU(2) ASL FIELD THEORY 

Here, we shall develop some of the details of the field 
theory of the SU(2) algebraic spin liquid in the honey- 
comb Hubbard model. We first work through the contin- 
uum limit, and write the Lagrangian in a form that man- 
ifestly exposes the Sp(4) symmetry discovered by Ran 
and Weni2£ Next, we enumerate the action of the micro- 
scopic symmetries on the fields of the low-energy theory. 
Finally, we discuss the 5-component vector of "compet- 
ing orders," in which the Neel vector and two-component 
VBS order parameter are unified together by the Sp(4) 
symmetry. 



1. Continuum limit and Sp(4) symmetry 
Consider the mean-field Hamiltonian 



(rr') 



(Bl) 



where the sum is over near-neighbor honeycomb bonds, 
and tp r , which carries the SU(2) gauge charge, is defined 
in Eq. (|4U|) . We will take the continuum limit and catalog 
the action of the microscopic symmetries on the contin- 
uum fields. 

As in Sec. [I] and Fig. [TJ we use the 2-site unit cell 
labeled by (R,i), with i = 1,2, and define the Fourier 
transform 



(B2) 



where N c is the number of unit cells in the lattice. The 
Hamiltonian then becomes 



f ipki 



where 



[H(k)] n = [H(k)} 22 = 0, 



(B3) 



(B4) 
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and 



[H(k)} 12 = [H(k)}* 21 = 1 



e 



— ik-a,2 



ik-(ai— 0,2) 



(B5) 



The Dirac nodes are at k = ±Q, where Q = (4w/3)x. 
Putting k = ±Q + q and expanding to first order in small 
q we find 



H(Q 



and 



\/3 , 1 o\ 

— (QxT +q y r ) 



H{-Q + q) = — (fer 1 - q y r 2 ). 



V>Q+q,2 



We define the continuum spinors 

01 (?) 

and 

02 (<?) ^ 



(B6) 



(B7) 



(B8) 



(B9) 



Note that (p a (with a = 1,2) is a 4-component object. 
As discussed in Sec. IVI Bl it is convenient to act on this 
space with tensor products of two different kinds of Pauli 
matrices. The t % Pauli matrices act on the Lorentz space, 
and the \x % Pauli matrices act in the SU(2) gauge space. 
The continuum Hamiltonian takes the form 



H c = v 



d 2 q 



(2nf 



(BIO) 

where v is the velocity. To make this look like the con- 
ventional Dirac Hamiltonian we define 



Ifl = T 3 (^l 
(fi 2 = T (pi, 



and we have 



d 2 q 



tfliqxT 1 + q y T 2 )ip a . 



(Bll) 
(B12) 



(B13) 



It is useful to write down again the relation between 
the continuum and lattice fields: 



Vi(<?) 



and 



4>-Q+q.2 
ll>-Q+q,l 



We also define the eight-component object 

<r» = [ ^ 

f2 



(B14) 



(B15) 



(B16) 



This introduces an additional SU(2) flavor space - we 
denote the Pauli matrices acting there by zA So, for 
example, 



v 1 ® = 



92 



(B17) 



We now write down the imaginary-time Lagrangian 
density 

Co = ^ a [d -i{T 1 d l + T 2 d 2 )\^ a (bis) 

= if a [ - ijfidfj] if a (B19) 
= $[-47^]*, (B20) 
where we have set v = 1 for convenience, and defined 

7m = (^WA (B21) 

and 



Va 



(B22) 
(B23) 



The coupling to the SU(2) gauge field a* is introduced 
by replacing <9 M with the covariant derivative 



(B24) 



resulting in the fermion kinetic energy part of the effec- 
tive Lagrangian density 



£ cS = $[- i^D = i p]$, 



(B25) 



where we have introduced the Feynman "slash" notation 

p w 

We now proceed to put C e s into a form that manifestly 
exposes the Sp(4) symmetry of Ran and Wen^ First, we 
define another kind of complex conjugate spinor 



$ = (in 2 ){iT 2 )<f> 7 



(B26) 



The form of <t takes advantage of the fact that the spin- 
1/2 representation of SU(2) is self-conjugate, so that $ 
and <E> transform identically under both Lorentz rota- 
tions and gauge transformations. For example, under 
a Lorentz rotation, <£> — * e IQAl7f '$, and also $ — > e lQ ^ 7fl $. 

This property suggests that we define yet another ob- 
ject, this one with sixteen components: 



T 



We also define a conjugate of T by 

T = T T (zr 2 )( iM 2 ). 



(B27) 



(B28) 



This has introduced yet another SU(2) space into our 
formulation; we use <, 1 Pauli matrices to act in this one. 
For example, 



(B29) 
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The advantage of introducing all this notation is that 
it allows us to write £ c g in a manifestly Sp(4)-invariant 
form. Some algebra verifies that, up to integration by 
parts, 



-\r^ FjT. 



c en = -^repr. (B30) 

This form is invariant under the global rotation 



MX, 



(B31) 



where M should be thought of as a 4 x 4 matrix acting in 
the tensor product of the v % and <; 1 Pauli matrix spaces. 
(So, for example, M commutes with // and r l Pauli ma- 
trices, but not, in general, v l and The condition that 
this transformation leaves C c g invariant is 



M T <?M 



(B32) 



This is precisely the statement that M is an Sp(4) matrix, 
so we have shown that £ c ff is invariant under global Sp(4) 
rotations. 

One point remains to be addressed. Because our repre- 
sentation of £ c ff is redundant to a certain extent, we have 
to check that the Sp(4) is a real symmetry and not just 
an artifact of the formalism. This can be done by rewrit- 
ing the Sp(4) Noether current in terms of the original 
fcrmion operators, and checking that it is non-vanishing 
and that its components are linearly independent. The 
Noether current is 



1. 



-T 7 ^ 2 T A T, 



(B33) 



where T A is one of the ten 4x4 matrix generators of 
Sp(4). Some straightforward algebra - or, alternatively, 
a group-theoretic analysis - shows that all components of 
J A are nonzero and linearly independent, and therefore 
the Sp(4) is a genuine global symmetry. 



2. Microscopic symmetries 

I now enumerate the action of the microscopic symme- 
tries on the continuum fields (in real space). The results 
quoted here can be derived starting from the definition 
of ip a in Eqs. (IB14IIB15|) . We need to consider the space 



group symmetries of the honeycomb lattice, spin-rotation 
symmetry and time-reversal. Pseudospin symmetry (in- 
cluding the particle-hole transformation) does not act on 
the spinon sector, so we do not discuss it here. 

First, we discuss the space group symmetries and time- 
reversal. (Spin rotations will be dealt with separately.) 
For each symmetry, we give its definition in terms of the 
microscopic spinon fields, then quote its action on the 
continuum fields T and T. We include enough space 
group operations to generate the entire space group. 

Translations. We consider translations by the basis 



vectors a\ and ct2. In general, Tr : ipR>.. 
the continuum we have 



R+R ,i ■ 



111 



T ai : T 
T ai : T 



T a2 :T 



T a2 :T 



exp(^V)T 



Texp (— i/V 
V 3 

'2m 



exp 



/2tt 
V~3 



Texp f^i/V 



(B34) 
(B35) 
(B36) 
(B37) 



Horizontal reflection. This operation is defined by TZ X 
■i/jRi -» ipRH, where R' = (-R x ,R y ). We have 



T(r) 
T(r) 



-cWT(r') 
T(r')[- 



?W1 



(B38) 
(B39) 



Vertical reflection. Here, we have !Z y : tpai ~> ^ 3r, pR'2 
and lZ y : %})R2 — > (X ipR'i, where R = {R x ,—R y ). The 
result in the continuum is 

K y : T(r) -» -? 3 i/ 3 rV 3 T(c') (B40) 
^:T(r) ^ T[-cVrV 3 ]. (B41) 

Rotation about plaquette center. We consider a coun- 
terclockwise rotation by 7r/3 about the plaquette center 
located at (1/2, 0). The action on the fermion fields is 

Rtt/3 ■ IpRl -> M 3 V'R'+a 1 -a 2 ,2 (B42) 
K/3 ■ 4>R2 -» ^IpR'l, (B43) 

where R' is the image of R under a rotation by 7r/3 about 
the origin. In the continuum, 



J 



R*/3 ■ T(r) 
K/3 ■■ T(r) 



<r 3 exp 



f (^V + ^)]exp(^yT(r') 



-T(r')«? 3 exp 



in ,1 [ 3 \/3 



exp i 



(B44) 
(B45) 



Time reversal. It is simplest to define time reversal in terms of the original /-fermions; we have T : f r 
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{i<y 2 ) a pf r p. It is important to remember that time rever- 
sal is an anti-unitary operation. Time reversal operates 
on the continuum fields as follows: 



T : T(r) 
T : T(r) 



-rWT(r) 



(B46) 
(B47) 



Now we consider spin rotations. As these make up part 
of the Sp(4) symmetry, we proceed somewhat differently 
from the symmetries above. First it is useful to provide 
some more detail about the Sp(4) symmetry itself. A 
general Sp(4) matrix can be expressed as M — e lX T , 
where A — 1, . . . , 10 labels the Sp(4) generators, and the 
condition M T <, 2 M = <; 2 implies 



2rpA 



(T ) q + fT 



0. 



We choose the following basis for the T A 
T 



l 2 2 2 2 J 



(B48) 



(B49) 



The conserved charges Q A = J d 2 r J A (r) generate Sp(4) 
rotations. Formally this means that 



„iX 



Q Te - 



T, 



(B50) 



where the objects in this equation should be interpreted 
as quantum operators, not as Grassmann variables. Ac- 
tually the form of the right-hand side of this equation 
is not completely obvious; depending on the normaliza- 
tion chosen for the Noether current, we could have a dif- 
ferent coefficient inside the exponential. To check that 
Eq. (|B50p is correct, it is enough to evaluate the left- 
hand side for some particular choice of X A ; this is just a 
matter of some algebra. 

To determine the action of spin rotations on T, the 
above discussion shows that we need only identify the 
spin density among the Sp(4) conserved densities J A . 
The correct identification is 



S x (r) 


— J o 




S y {r) 


- T A 


rpA 


S'(r) 


- T A 


rpA 



"2 ? 



(B51) 
(B52) 
(B53) 



So, denoting by i? s (a) a spin rotation by an angle |a| 
about the o:/|q:| axis, we have by Eq. (|B50|) 



R a (6x) : T 
R s (9y) : T 
Rs(0z) : T 



exp 



T 



(B54) 
(B55) 
(B56) 



3. Unification of competing orders 

Here, we construct the gauge-invariant, Lorentz-singlet 
fermion bilinears in the SU(2) ASL. There are six such 
bilinears; five of them transform as the SO (5) vector irre- 
ducible representation of Sp(4), and the other is an Sp(4) 
singlet. Ran and Wen2& found that the correlations of the 
SO(5) vector bilinears are enhanced by gauge fluctuations 
(see Sec. IVI B[) . We show here that the SO(5) vector is 
a unification of the Neel vector, and the two-component 
order parameter for the VBS state pictured in Fig. [4] We 
also define the VBS order parameter microscopically in 
terms of spin operators. 

A general gauge-invariant, Lorentz-singlet fermion bi- 
linear has the form 



B 



w 



TWT, 



(B57) 



where W is a 4 x 4 matrix acting in the tensor product 
space of the ^ and v l Pauli matrices. We require Bw to 
be Hermitian, which puts the following condition on W: 



W 



(B58) 



There are 16 linearly independent matrices satisfying 
this condition. We wish to decompose the resulting 16 
bilinears into irreducible representations of Sp(4). To 
do this, note that under an Sp(4) rotation T — > MY, 
W -> M T WM. This transformation preserves the sym- 
metry of W, so we can consider symmetric and antisym- 
metric W matrices separately. Furthermore, tv{q 2 W) is 
invariant under Sp(4) rotations, so it must be constant 
within each irreducible representation. We are therefore 
led to the following bases of W matrices for the irre- 
ducible representations 



wis 



{^VV^kV^V} (B59) 

{C 2 } (B60) 
{<; V , 1 , j/ , z/ , k 1 , k 1 j/ 1 , k V , 

K 3 ,kV,K 3 i/ 3 }. (B61) 



The first two of these consist of antisymmetric W matri- 
ces; W\ s is the SO(5) vector representation, and W\ s 
is an Sp(4) singlet. It can be seen that all bilinears Bw 
with the symmetric W matrices of Ws vanish identically; 
this is a consequence of Fermi statistics and the redun- 
dancy of our representation for the fermion fields. This 
can be seen by showing, for example, that zT^T = 0; all 
the other bilinears obtained from Ws must then vanish 
by Sp(4) symmetry. 

Therefore we are left with only six bilinears. We orga- 
nize these by defining 

N A = (t^ 2 (-^V 2 )T,T<; 2 ( 1 ; 2 ^ 2 )T,T 1 ; 2 (i, 3 )T, 

T<; 2 (^)T,TV 2 (<;V)T), (B62) 

where A = 1, . . . , 5, and 

fi = T<r 2 T. (B63) 
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N A is the SO (5) vector and fl is the Sp(4) singlet. 
Next, we define 



n = (V 1 , V 2 , V 3 ) 



and 



(N 4 ,N 5 



(B64) 



(B65) 



It is a straightforward exercise to use the results of the 
previous section to show that n transforms like the Neel 
vector under all microscopic symmetries, and v trans- 
forms like the two-component order parameter for the 
VBS state pictured in Fig. [U 

We now define the VBS order parameter microscopi- 
cally in terms of spin operators. It is analogous to the 
order parameter used by Read and Sachdev to charac- 
terize VBS states on the square latticed We define the 
order parameter as 



V"VBS 



E ex p( 



-R ai 



R 

e &ri ■ '-'H+ai-02,2 

„ — 2-Ki/3 q G 



The relation between this object and v is 



J d 2 x(v 1 (x) + iv 2 {x)) ~ Wbs- 



(B66) 



(B67) 



Now consider the pattern of VBS order show in Fig. [4] 
and suppose that, on the strong bonds, (S r ■ S r <) = c%, 
while on the weak bonds, (S r ■ S r >) — ci- Then it can be 
shown that, in this state, 



1 

iV 



(V'- 



3(ci - c 2 ) 



(B68) 



where N c is the number of unit cells in the crystal, which 
should of course be taken to infinity to obtain the ther- 
modynamic limit. The phase 9 can take on the three 
values 9 — 0,±27r/3, which correspond to the three dis- 
tinct VBS ground states. 

As an aside, it is interesting to note that, if the roles of 
strong and weak bonds are reversed, the "bond-ordered" 
state of Fig. 0] becomes a plaquette state, and (V'VBS ) 
simply changes sign. Therefore, the pattern of symmetry 
breaking is identical in the bond-ordered and plaquette 
states, and it is not at all obvious that these two states are 
really sharply distinct phases. On the square lattice, the 
Read-Sachdev VBS order parameter also contains both 
the columnar and plaquette VBS states, but these states 
do not have the property of being related by interchang- 
ing strong and weak bonds, and are clearly characterized 
by different patterns of symmetry breaking. We note that 
a study of the honeycomb lattice quantum dimer model^ 2 - 
found a first-order transition between bond-ordered and 
plaquette states, which does not exclude the possibility 
that these states are in fact different limits of the same 
phase. 



APPENDIX C: FORMAL EQUIVALENCE OF 
SU(2) AND U(l) MEAN-FIELD CALCULATIONS 

Here, we demonstrate the formal equivalence of the 
mean-field calculations for the SU(2) and U(l) slave ro- 
tor formulations of the Hubbard model. This equivalence 
holds for the restricted class of mean-field ansatz for the 
SU(2) formulation given in Eqs. (|91fl93|) . We may there- 
fore use the results of Ref.[l6|for the honeycomb Hubbard 
model in our discussion of the mean-field phase diagram 
in Sec. IVBl 

Below, we shall consider saddle points with a T = 0. In 
that case, the ground state energy in mean-field theory 
is given by 



771MFT 

(r,r') 



t r A\x F rr>\\X Z rMUyU r z r ) 

+ Eg + Eg, (CI) 



where Eq and Eq are the ground state energies of the 
fcrmions and bosons, respectively. The fermion energy 
can easily be obtained by taking the expectation value of 
the fermion mean-field Hamiltonian: 

Eg = E t rr ,\ X U^{F r ,U^Ft) (C2) 

(r,r') 

= - E \trr'\\Xrr'\\Xrr'\tr(U r F r U r /)- (C3) 
(r ,r') 

This cancels the first term of Eq. (|CJ1[) , and 



E, 



MFT 



= - e -v 



(C4) 



To demonstrate the equivalence of the present mean- 
field calculation and that of Ref. [l6|, it is first useful to 
plug our ansatz into the saddle point action. The result 
for the rotor boson term is 



~ Sz = I dT [T, z U-^ T + ~x r 

E \U r i\\Xrr>\{^'zUZrH 
(r,r') 



— d T + £7 I 



(C5) 
(C6) 



where in the second equality we have transformed to a 
basis labeled by 7 that diagonalizes the quadratic form, 
with eigenvalues e 7 . Explicitly, 



— ^ t f~/{ r )z~fi- 



(C7) 



The fermion part of the saddle-point action is 

Sf — dr I" E frgd-rfra 
J r 

E ^'|X^l(e^7ra/wa+H.C.)].(C8) 
(r,r') 
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Using Eq. (|C4[) , the mean-field ground state energy is 



saddle-point equations have the form 



E, 



MFT 



The saddle-point equations become 

1/2 = (|z rl | 2 ) = (| Zr2 | 2 ) 
= (ipt<nP r ) 



\Xr 



2(z* 1 Zr'l) 



2{z* 2 z r , 2 ) 



(C9) 



(CIO) 
(Cll) 
(C12) 
(C13) 



and the remaining equation, lm(tr(d T Zj,crZ r )) = 0, is 
automatically satisfied by saddle points of the form con- 
sidered. 

Now we will briefly introduce the U(l) formulation, 
and then show that it reproduces the same saddle-point 
equations and ground state energy. Following Ref.[l6|, the 
partition function of the Hubbard model can be written 

Z = J VfVfVXVX*VhVXexp(- Su(i)), (C14) 
where <SYj(i) — / ^ T ^u(i) ( T ) an d 

LV(1) = 7^77 ^ \ ( id r + h r )X r \ 2 + y^Jra(d T + ih r ) f ra 
r r 

- iJ2h r + iJ2x r (x;x r -i) 

r r 

— (t rr i f ra f r > a X*X r i + H.c.) . (C15) 

(r,r ; ) 

Here it is important to note that we have taken the 
coefficient of the on-site repulsive interaction to be U' 
rather than U. It will turn out that the U(l) and 
SU(2) mean-field calculations are equivalent provided 
that U' = (4/3)17. 

The hopping term is decoupled using a pair of complex 
fields rj rr i and rj r > r , resulting in the following contribu- 
tions to the action: 



5 u(i) 



S. 



U(l) 

tx 



dr \trr'\(\V 

— \t rr '\{X^.rj rr > X 

(r,r') 
X r fT} r f r X r ) 



it' I \Vr'r | ) 



s, 



U(l) 

tf 



J ^ ^ t rr ' {^fraVrr' fr'ct 



(r,r' ) 
fr'aVr'rfrot) 



(C16) 



(C17) 



(C18) 



We can now look for saddle points by varying the ac- 
tion with respect to r] rr , ij r > r , X r and h r . The resulting 



rj* r , — (X*X r i) 

Tj rr i sign(^ rr ' ) {fra fr'i 

i = (x;x r ) 



(C19) 
(C20) 
(C21) 



= —hr + i{f ra fra) + ^((d T X*)X r ),(C22) 

where we have not explicitly written the equations in- 
volving rj r i r , which are similar in form to those involving 
rj rr i . We now write 



Vrr' \Xrr'\^ rr & 



7] T ' r — IXrr'l^ rr G 

and make the analytic continuation 



(C23) 
(C24) 

(C25) 





-» ia+ r , 


(C26) 


X r 




(C27) 


h r 


— » ih r + 5h r . 


(C28) 


We look for saddle points where h r = 0. 


Then, defining 


1 Xt")"' 1 


EE |Xrr'|e- a ^ 


(C29) 


1 X-rr' 1 


EE IXrr'le 5 -' 


(C30) 


~ X 

CLj.y.1 


— UJ'P'p! | (Xj"p' 


(C31) 


u rr' 




(C32) 



the saddle-point equations are 

\Xrr-\e 1 ^' = sign(t rr ,)(f ra f r , a ) (C33) 
\xL\e' ra ^' = {X* r X r ,) (C34) 



1 = (X:X r ) 

1 {fra fra) 

= lm{(d T X*)X r ). 



(C35) 
(C36) 
(C37) 



As in the SU(2) case, Eq. (|C37|) is trivially satisfied. The 
expectation values are to be calculated using the saddle- 
point action <Stj(i) — / ^ T -7u(i) ( T )i where 



-U(l) 



2 J2 1^' 1 I \Xrr' I COs(a£, - ) - A, 

(r,r') r 

+ frMra + £ K ( ~ Tj^r + X r ^X r 

r r 

- J2 trr' \xlr> I (e _<3 ^ /m/r-'a + H.C.) 
(r,r ; ) 

£ | W | | x £, | {X*e<^ Xr , + c.c.) . (C38) 

(r,r') 

If we identify \x f rr ,\ = \Xrr'\ and «f r ' = ^r" the 
fermion part of Eq. (|C38|) is identical to Sf of the SU(2) 
formulation. All expectation values of fermion operators 
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are then equal in the two formulations, and by Eq. (IC33|) 

we Should also identify \Xrr'\ = \Xrr'\ an ^ 0>tt' = ^rr>; 

So far we have shown that the saddle-point Eqs. (|C33|) 
and (|C36p are identical to their counterparts in the SU(2) 
formulation. 

The part of Eq. (|C38|) involving X r can be written 



L 



u(i) _ 
x ~ 



c 7 ) 7 ' 



where 



(C39) 



(C40) 



and (|C35p . identical to their SU(2) counterparts. To see 
this, we evaluate the equal-time correlator (in the SU(2) 
mean-field theory) 

(z* rl z rll ) = ^^/ 7 *(r)/y(r')(^^yi) (C43) 

7 7' 

, 2n Z.( 3 /2C/K + e7 - ^ > 

Similarly, in the U(l) formulation, putting U' — (4/3)U, 
we have 



and / 7 (t*) and e 7 are the same eigenfunctions and eigen- 
values as in Eqs. (|C6|) and (|C7[) . Proceeding as above, 
the ground state energy is 



E, 



MFT 



'2U> 



7 



(C41) 



This can be made equal to the ground state energy cal- 
culated in the SU(2) formulation, Eq. (|C9|) . by choosing 



(C45) 



_ wo fC46) 



where in the last line we have changed variables by ui v 
2tv n . Therefore we have shown that 



U' 



:U. 



(C42) 



{X*X r ,) = 2{z* rl z r/1 ) = 2{z* r2 z r , 2 ), 



(C47) 



We need to check that this choice for U' will also 
make the remaining saddle-point equations, Eqs. (|C34j) 



which implies that the remaining saddle-point equations 
are identical in the two formulations. 
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